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SECTION  1 
INTRODUCTION 

.  NASTRAN  is  a  large,  general  purpose,  finite  element  computer  program 
used  primarily  for  linear  structural  analysis.  In  its  ten  years  of  existence, 
NASTRAN  has  achieved  the  distinction  of  being  the  program  with  the  most  users 
worldwide,  most  documentation,  most  user  conferences,  most  cited,  and  probably 
the  most  expensive  to  develop. 

The  four  objectives  of  this  study  were  to: 

.Assess  the  current  status  (usage,  documentation,  trends)  of  public- 
domain  versions  of  NASTRAN,' 

l  / 

^Evaluate  its  program  architecture,  and  compare  it  with  some  other 
general  purpose  codes^ 

, /purvey  its  functional  description,  and  comment  on  its  capabilities/ 
limitations  and  such  items  as  matrix  operations,  eigenvalue  extraction 
schemes,  element  library,  material  properties,  etc/ 

.Perform  advanced  evaluation  exercises  on  some  selected  structural  elements, 
and  test  their  convergence,  completeness,  and  efficiency. 

Obviously,  it  is  impossible  in  a  study  of  this  limited  scope  to  evaluate 
comprehensively  all  aspects  of  NASTRAN.  Therefore,  at  the  outset,  it  was  decided 
to  focus  attention  at  the  following  areas  in  the  first  half  of  the  year-long 
study:  user  survey,  documentation  review,  program  architecture,  functional  descrip¬ 
tion,  and  selected  verification  exercises.  The  second  half  of  the  study  concen¬ 
trated  on  the  advanced  evaluation  exercises,  designed  to  assess  element  quality 
and  limitations  (see  Section  7). 

There  are  literally  thousands  of  NASTRAN  users.  For  convenience  in  this 
study,  we  shall  divide  them  into  four  categories: 


_ CATEGORY _  _ DESCRIPTION _ 

1.  Expert  This  person  knows  NASTRAN  inside  out.  He  is  familiar 

with  all  the  documentation,  including  the  NASTRAN 
Programmer 's  Manual.  He  is  adroit  with  the  Direct 
Matrix  Abstraction  Program  (DMAP) .  He  has  used 
NASTRAN  in  a  variety  of  problems  for  probably  close 
to  ten  years,  and  is  recognized  in  his  organization 
as  the  NASTRAN  expert,  the  person  to  whom  everyone 
comes  for  help. 


CATEGORY 


DESCRIPTION 


2.  Experienced  This  person  knows  NASTRAN  well,  and  has  probably  used 

the  code  for  five  to  seven  years.  He  is  competent  in 
finite  element  theory  and  modeling,  and  most  likely, 
has  used  several  other  codes.  He  knows  how  to  plot 
models  to  his  liking  with  minimum  effort.  He  has  used 
several  of  the  NASTRAN  rigid  formats  and  DMAP,  and  knows 
the  input/output  well.  This  person  is  usually  a  senior 
engineer  in  a  large  company. 

3.  Common  User  The  common  NASTRAN  user  has  generally  used  the  code  for 

two  to  five  years  and  is  familiar  with  rudimentary  finite 
element  theory.  He  is  likely  to  have  used  one  rigid 
format  (e.g.  statics)  almost  exclusively,  to  analyze 
different  designs  of  one  product.  This  person  finds  the 
concept  of  DMAP  very  difficult. 

4.  Novice  The  new  NASTRAN  user  has  zero  to  two  years  experience.  He 

is  becoming  familiar  with  NASTRAN  modeling,  input,  and 
output.  He  needs  guidance  in  modeling  techniques,  works 
for  the  experienced  or  common  user,  and  is  apt  to  make 
_ mistakes  in  element  connectivity  and  boundary  conditions. _ 

Throughout  this  report,  certain  comments  and  conclusions  shall  be  addressed  to 
each  category  of  users.  Most  NASTRAN  users  fall  in  Categories  2  and  3. 

Unlike  most  current  general  purpose  finite  element  computer  programs  de¬ 
veloped  by  individuals,  companies,  and  universities,  NASTRAN  was  originally 
conceived  by  and  developed  under  NASA's  sponsorship  in  1964-1969.  After  con¬ 
sulting  with  an  ad-hoc  committee  of  finite  element  experts  from  the  aerospace 
industry,  NASA  established  a  NASTRAN  project  management  office  which: 
administered  the  software  development  contracts,  released  and  updated  periodically 
the  voluminous  documentation,  established  a  systematic  error  reporting  and 
correction  procedure,  and  shaped  the  contents  of  the  code.  This  responsibility 
was  shifted  in  1979  to  the  Computer  Software  Management  and  Information  Center 
(COSMIC)  at  the  University  of  Georgia.  NASA  was  therefore  the  "code  developer" 
of  NASTRAN.  This  short  explanation  on  NASTRAN  development  history  serves  to 
point  out  the  fact  that,  unlike  other  competing  codes,  COSMIC/NASTRAN  currently 
does  not  offer  a  "hot-line  service"  to  answer  promptly  user  questions  on  the 
telephone  (see  end  of  Section  2) . 

The  future  success  and  usefulness  of  NASTRAN  (as  well  as  any  other  structural 
analysis  software)  is  intimately  related  to  the  current  and  future  developments 
in  computers.  For  some  excellent  recent  reviews  and  projections  on  this  subject. 


the  reader  is  referred  to  References  1  to  6.  These  articles  review  computer 
developments  in  the  past  thirty  years  and  project  trends  into  the  1980's. 

Topics  which  are  covered,  most  of  which  will  have  an  effect  on  the  future  of 
NASTRAN,  are:  operating  systems,  high-level  languages,  time  sharing,  CAD/CAM, 
minicomputers,  interactive  terminals,  new  memory  chips,  virtual  and  bubble 
memories,  disc  files,  distributed  processing,  and  supercomputers  ("number 
crunchers").  The  computer  industry  is  on  the  threshold  of  very  large  scale 
integration  (VLSI)  chip  technology  and  bubble  memories,  both  of  which  appeared 
on  a  commercial  scale  in  1979.  The  future  impact  of  these  developments  on 
structural  analysis  software  is  difficult  to  predict.  In  view  of  the  current 
proliferation  of  minicomputers  on  the  structural  analysis  scene,  it  was 
decided  at  the  end  of  the  NASTRAN  advanced  evaluation  (Section  7)  to  comment 
briefly  on  minicomputer  versions  of  NASTRAN  and  also  on  pre-  and  post-processors 
now  commercially  available  for  NASTRAN.  The  real  worth  of  a  code  lies  in  its 
maintenance  and  support. 
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USER  SURVEY  RESULTS 

This  section  summarizes  the  findings  of  a  telephone  survey  of  approxi¬ 
mately  25  NASTRAN  users,  most  of  whom  fall  in  Categories  2  and  3  and  work  at 
large  aerospace  companies.  This  survey  was  by  no  means  a  comprehensive  picture 
of  NASTRAN  usage;  rather,  it  was  basically  a  regional  consensus  which  reflected 
each  user's  own  opinions  and  experiences.  The  aims  of  the  survey  were  to  obtain 
a  quick  feel  of  current  NASTRAN  usage,  obtain  user  comments  and  experiences,  and 
solicit  suggestions  for  improvement. 

The  main  problem  with  the  survey  was  to  find  current  users  of  public- 
domain  versions  of  NASTRAN.  In  fact,  only  a  few  users  of  the  latest  release 
(Level  17.5)  of  COSMIC/NASTRAN  were  found.  Many  users  had  learned  NASTRAN  using 
the  COSMIC  version,  but  had  since  switched  to  commercial  versions  (such  as 
Universal  Analytics  Inc.'s  UAI/NASTRAN  and  MacNeal-Schwendler  Corporation's 
MSC/NASTRAN) .  In  the  past,  many  government  agencies  in  awarding  contracts  had 
dictated  that  the  vendor  use  NASTRAN  for  structural  analysis,  without  specifying 
the  particular  public  or  commercial  version.  Therefore,  the  vendor  could  choose 
the  one  he  preferred.  The  overwhelming  reason  for  switching  to  commercial  ver¬ 
sions  like  UAI/NASTRAN  and  MSC/NASTRAN  was  better  user  support  for  day-to-day 
problems.  Other  important  reasons  were  significant  cost  improvements  and  more 
efficient  solution  algorithms. 

The  following  statements  summarize  some  major  findings  of  the  survey.  Several 
interesting  comments  are  amplified  later. 

.There  still  exists  a  large  group  of  goverment  and  industry  users  who  prefer 
the  public-domain  versions  of  NASTRAN,  and  who  like  to  modify  the  code  to 
suit  their  own  needs. 

.The  biggest  single  complaint  of  COSMIC/NASTRAN  users  is  the  lack  of  a  hot¬ 
line  service  for  prompt  support.  (See  comment  at  the  end  of  this  section) . 

.NASTRAN  users  generally  think  of  the  code  favorably,  and  consider  their 
NASTRAN  experience  as  a  valuable  asset  in  their  career  development. 

.Most  users  confine  their  analyses  to  rigid  formats  and  perform  linear  elastic 
analysis.  A  few  of  the  sophisticated  users  create  their  own  DMAP  sequences, 
but  this  appears  to  be  quite  rare. 

.Another  big  complaint  is  cost.  NASTRAN  is  rated  too  costly  for  small  and 
medium-sized  jobs.  Efficiency  and  run-time  improvements  are  reported  in 
several  sources  (for  example,  References  14  and  40)  comparing  MSC/NASTRAN 
and  UAI/NASTRAN  to  COSMIC/NASTRAN,  but  the  present  study  will  not  evaluate 
this  aspect. 
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.NASTRAN  documentation  is  bulky,  and  difficult  to  read  and  use  (see 
Section  3) . 

.There  is  a  need  for  an  internal  mesh  generation  capability. 

•Many  users  wished  NASTRAN  had  nonlinear  material  and  geometry  capabilities. 
(Its  present  nonlinear  capability  is  limited,  and  consequently  will  not 
be  evaluated  in  this  study) . 

Most  users  agreed  that  the  executive  system- functional  modules  approach 
in  NASTRAN  program  architecture  makes  for  a  very  flexible  and  powerful  program 
(see  Section  4).  They  use  a  few  rigid  formats  heavily,  and  rarely  have  need 
for  creating  DMAP  sequences.  Usually,  DMAP  alters  are  used  for  output  variations. 
The  program  can  be  used  for  matrix  manipulations,  again  a  feature  due  to  DMAP  and 
the  modular  concept.  This  feature  brought  out  another  interesting  observation 
on  the  program  architecture.  NASTRAN  was  conceived  and  developed  ten  to  fifteen 
years  ago,  when  computer  core  was  expensive  and  the  modules  were  created  to  limit 
the  dependence  on  core  size.  Today,  core  is  plentiful  and  cheap  (and  rapidly 
becoming  cheaper),  virtual  memories  are  commonplace  and  bubble  memories  are 
emerging.  Most  contemporary  programs  are  manipulated  by  the  machine  operating 
system  rather  than  an  executive  system  internal  to  the  program  (such  as  NASTRAN) . 
This  NASTRAN  program  architecture  prevents  the  user  from  communicating  inter¬ 
actively  with  the  program. 

In  the  area  of  error  and  bug  correction,  the  user's  views  depend  largely 
on  the  size  of  his  organization  and  whether  there  exists  a  competent  in-house 
staff  to  correct  bugs  and  implement  new  elements  or  modules.  The  small  or 
medium-sized  user  finds  it  virtually  impossible  to  get  things  corrected  on  his 
own  on  COSMIC/NASTRAN ,  and  will  very  likely  opt  for  a  commercial  version.  The 
.  ger  user  has  the  computing  and  manpower  resources  necessary  for  corrections 
and  implementations,  and  may  prefer  COSMIC/NASTRAN  for  its  flexibility.  Appar¬ 
ently,  the  released  versions  of  the  program  from  COSMIC  have  not  been  thoroughly 
tested  and  verified,  causing  many  reported  system-type  problems  and  errors.  When 
a  user  discovers  a  bug,  his  only  course  of  action  (other  than  correcting  it 
himself)  is  to  report  it  to  the  current  maintenance  contractor.  Then,  the  govern¬ 
ment  decides  on  the  priority  of  correcting  various  reported  errors,  and  the 
maintenance  contractor  may  not  be  able  to  correct  the  error  until  the  next  released 
level.  Another  problem  is  the  development  effort  required  for  the  different 
machine  versions  of  a  given  NASTRAN  release  level.  A  level  is  usually  developed 
on  one  machine,  and  the  testing  reliability  and  quality  assurance  may  not  be  as 
high  on  another  machine.  * 
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Users  complained  about  possible  errors  in  the  capabilities  dealing  with: 
substructuring,  cyclic  symmetry,  differential  stiffness  with  thermal  loading, 
axisymmetr ical  analysis,  piecewise  linear  analysis,  eigenvalue  extraction  methods 
and  their  accuracy  in  computing  rigid-body  modes  for  large  problems,  and  the 
QUAD 2  element.  Capabilities  and  options  not  presently  in  NASTRAN  but  desired 
by  many  users  are:  incompressible  elements,  composite  materials,  strain  recovery, 
mesh  generation,  timing  estimates,  non-zero  initial  conditions  for  modal 
analysis,  alternate  time  integration  schemes,  and  nonlinear  geometric  and  material 
analyses. 

Commentary 

A  very  interesting  comment  on  NASTRAN  maintenance  was  made  by  several  users. 
COSMIC/NASTRAN  was  first  released  to  the  public  in  1969,  and  has  been  updated 
through  continuous  government  contracts.  Each  new  updating  contract  is  bid  on 
by  private  companies.  Then,  the  successful  bidder's  responsibility  is  to  define, 
initiate,  integrate,  and  quality-assure  all  work  done  on  the  program.  Therefore, 
each  new  level  released  to  the  public  is  a  direct  reflection  of  the  performance  of 
the  maintenance  contractor  and  subcontractors.  In  the  past,  continuity  and 
coordination  of  efforts  on  some  occasions  could  have  been  better. 

Currently,  user  support  on  COSMIC/NASTRAN  is  coordinated  by  Dr.  Robert  L. 

Brugh,  NASTRAN  Project  Manager  at  COSMIC,  University  of  Georgia  at  Athens, 

Telephone:  (404)  542-3265.  He  is  responsible  for  providing  answers  to  user’s 
questions  or  problems,  and  coordinates  closely  with  Computer  Sciences  Corporation, 
the  current  COSMIC/NASTRAN  maintenance  contractor.  COSMIC  maintains  and  updates 
the  NASTRAN  Software  Problem  Report  (SPR)  log,  which  can  be  requested  by  any  user. 
Several  recent  reports  from  COSMIC/NASTRAN  users  with  specific  questions  and 
problems  (and  who  decided  to  use  this  channel  of  communication  with  Dr.  Brugh 
and  COSMIC)  have  been  favorable.  COSMIC/NASTRAN  users  are  therefore  urged  to  contact 
COSMIC  promptly  when  they  discover  system  problems,  element  anomalies,  and  other 
bugs. 
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SECTION  3 

NASTRAN  DOCUMENTATION 

Some  general  comments  are  made  here  about  NASTRAN  documentation,  which 


primarily  consists  of-' 

.The  NASTRAN  Theoretical  Manual  (Ref.  7) 

.The  NASTRAN  User's  Manual  (Ref.  8) 

.The  NASTRAN  Programmer's  Manual  (Ref.  9) 

.The  NASTRAN  Demonstration  Manual  (Ref.  10) 

.NASTRAN  User's  Guide  (Level  17.5)  (Ref.  11) 


The  first  four  manuals  were  designed  for  reference  purposes,  and  are  bulky 
and  sometimes  difficult  to  use.  In  general,  the  documentation  in  these  four 
manuals  is  fairly  complete  though  of  uneven  quality.  Of  the  four  user 
categories,  the  expert  user  is  probably  the  only  one  to  have  examined  all 
these  thick  volumes.  It  is  often  difficult  to  find  information.  Examples 
of  this  difficulty  are  given  below: 


MANUAL _ "DIFFICULT-TO-FIND"  INFORMATION  (OR  MISSING) _ 

Number  of  nodes  and  degrees  of  freedom  for  each  element. 

User  hints  and  limitations  of  each  element. 

Section  5.8  (80  pages  long)  adequate  description  of  15  membrane 
and  bending  2-D  surface  elements,  but  element  names  are  missing 
in  half  of  the  subsections. 

Element  summary  table  showing:  classification,  name,  brief 
description,  nodes,  degrees  of  freedom,  material  options,  and 
element  developer  and  date. 

Time  and  cost  estimates  for  different  rigid  formats  and  elements. 

Tables  showing:  geometry  cards?  heat  transfer,  fluid,  mass, 
rigid  and  dummy  elements?  material  properties?  constraints  and 
partitioning?  loads  (static,  dynamic,  heat  transfer)?  problem 
control. 

Number  of  overlay  levels  and  system  flow  charts  for  the  IBM, 
UNIVAC,  and  CDC  computer  systems. 

Simple  problems  explaining  input  stream  line-by-line  to  the 
novice  user. _ _ _ _ 


A 


Theoretical 

Theoretical 

Theoretical 

User '  s 

User's 

User's 

Programmer ' s 
Demonstration 


Most  experienced  and  common  users  have  undoubtedly  read  only  those  sections 
in  the  manuals  necessary  for  their  day-to-day  use.  The  last  two  sections  in 
the  User’s  Manual  are  particularly  helpful  and  commendable  (since  they  do  not 
usually  appear  in  most  other  code  manuals):  the  well-catalogued  diagnostic 
messages  and  the  dictionary  of  NASTRAN  terms.  An  attempt  to  improve  many  of 
the  above-mentioned  difficulties  and  missing  items  in  the  four  manuals  of 
NASTRAN  documentation  is  made  in  the  NASTRAN  User's  Guide  (Level  17.5). 

The  User's  Guide  is  indispensable  to  all  users  of  COSMIC/NASTRAN .  It 
gives  a  concise  overview  of:  the  documentation;  NASTRAN  modeling;  descriptions 
of  the  various  rigid  formats;  program  architecture;  input  deck  description; 
executive  control,  case  control,-  and  substructure  control  decks;  printing  and 
plotting  features;  NASTRAN  implementation  and  input  decks  for  the  IBM  360/370 
(OS),  UNIVAC  1108/1110  (EXEC  8),  CDC  6000/CYBER  (NOS  and  NOS/BE)  computer  systems; 
NASTRAN  information  sources;  estimate  of  resource  and  time  requirements;  DMAP; 
and  application  examples  for  all  major  types  of  analyses.  This  618-page  guide 
is  an  admirable  attempt  to  place  in  one  document  most  of  the  useful  information 
a  user  will  ever  need.  The  following  items  are  especially  valuable:  a  table 
summarizing  NASTRAN  finite  element  names  versus  their  characteristics  (nodes, 
mass,  load  types  allowed,  heat  transfer  capability,  etc.);  a  general  discussion 
of  general  modeling  restrictions  in  NASTRAN;  a  summary  table  of  analysis  options 
for  dynamic  rigid  formats;  a  comprehensive  summary  table  of  modeling  options 
versus  each  of  the  20  rigid  formats  (15  displacement,  3  heat  transfer,  2  aero- 
elasticity)  ;  many  examples  of  executive  and  case  control  deck  setups;  a  13-page 
summary  table  of  bulk  data  options  (badly  needed  in  front  of  the  User's  Manual); 
a  clear  explanation  of  the  NASTRAN  overlay  system;  practical  guidelines  to 
estimate  cost,  rigid  format  resource  requirements,  CPU  time,  and  other  hardware 
requirements;  and  finally,  many  detailed  line-by-line  input  deck  explanations 
of  application  examples  illustrating  the  proper  use  of  all  major  rigid- format 
analyses.  Therefore,  the  User's  Guide  is  strongly  recommended  as  the  first 
NASTRAN  document  to  read  for  the  novice  user,  and  as  a  handy  reference  for  the 
common,  experienced,  and  expert  users  alike. 

Other  NASTRAN-related  documentation  which  may  be  of  interest: 

.NASTRAN  User's  Experiences  and  Colloquia  (First  through  Eighth,  1971- 
1979,  available  from  National  Technical  Information  Service,  Springfield, 
Virginia  22161) . 


.NASTRAN  Software  Problem  Report  (SPR)  Log  (periodically  updated  and  published 
by  COSMIC) . 

.NASTRAN  Newsletter  (free  subscription,  from  COSMIC). 

.Documentation  Error  Reports  (DER)  (COSMIC) . 

.NASTRAN  -  NASA  Structural  Analysis  (1979)  (COSMIC) . 

. MSC/NASTRAN  Basic  Training  Manual  (Ref.  12). 

. MSC/NASTRAN  Primer:  Static  and  Normal  Modes  Analysis  (Ref.  13). 

•A  Brief  Description  of  MSC/NASTRAN  (Ref.  14). 

f compares  MSC/NASTRAN  and  Level  16  COSMIC/NASTRAN ,  as  of  April,  1978; 

This  document  is  scheduled  to  be  updated  in  Spring  1981  and  will  compare 
MSC/NASTRAN  to  Level  17.5  Cosmic  NASTRAN 

.MSC/NASTRAN  Applications  Manual  (Ref.  53) 

.MSC/NASTRAN  Demonstration  Manual 

.MSC/NASTRAN  Handbook  for  Linear/Static  Analysis  (June  1980) 


SECTION  4 

PROGRAM  ARCHITECTURE 

The  key  concept  behind  NASTRAN  program  architecture  is  modularity.  When 
NASTRAN  was  originally  designed,  certain  design  criteria  were  stipulated 
(References  7,  9,  11),  such  as:  simplicity  of  input,  minimization  of  chances 
for  human  error  and  the  need  for  human  intervention  during  program  execution; 
functional  independence  of  solution  modules;  ease  of  program  modification  and 
extension;  versatility  and  adaptability  for  various  computer  systems;  restart 
capability,  etc.  The  resulting  NASTRAN  program  architecture  is  shown  in  Figure 
4-1,  which  illustrates  the  major  functional  subdivisions  and  their  interrelation¬ 
ships. 

The  two  major  components  are  the  Executive  System  and  the  Functional  Modules. 

The  Executive  System  is  the  heart  of  the  system.  It  executes  the  program  in 
two  phases:  the  preface  (where  a  basic  setup  is  performed) ;  and  the  execution  of 
the  DMAP  program  in  which  the  modules  are  controlled.  The  preface  does  extensive 
pre-processing  and  sorting  of  bulk  data  (this  function  of  NASTRAN  is  superior  to 
most  general  purpose  codes),  generates  initial  file  allocation  tables,  and  initial¬ 
izes  the  problem.  The  Executive  System  establishes,  protects,  and  communicates 
values  of  parameters  for  each  module.  It  allocates  system  files  to  all  data  blocks 
generated  during  program  execution;  a  file  is  "allocated"  to  a  data  block,  and  a 
data  block  is  "assigned"  to  a  file.  It  also  maintains  a  full  restart  capability 
for  restoring  program  execution  after  either  a  scheduled  or  unscheduled  interruption. 
The  general  philosophy  in  programming  Executive  System  routines  is  that  reliability 
and  efficiency  are  paramount  concerns.  Few  changes  were  anticipated  in  the 
Executive  System  as  NASTRAN  grew.  Therefore,  a  few  general  rules  were  imposed  on 
the  Executive  System  code,  resulting  in  a  sophisticated,  but  difficult  to  modify, 
system  (Ref.  11).  The  Executive  System  is  dependent  on  a  particular  computer 
and  its  operating  system.  It  comprises  approximately  10  percent  of  the  total  code. 

The  Functional  Modules  are  self-contained  subprograms  which  may  not  call, 
or  be  called  by  other  modules.  They  may  be  entered  only  from  the  Executive  System. 
Each  module  has:  data  blocks  (an  important  type  is  a  matrix),  subroutines  (which 
can  communicate  with  each  other  within  a  module);  parameters;  and  drivers.  A 
change  in  one  module  affects  no  others,  as  long  as  the  interface  with  the  Executive 
System  is  preserved.  No  module  can  directly  specify  or  allocate  physical  files. 


Figure  4-1  hastrah  program  architecture  and  functional  subvivisions 
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A  module  can  vary  in  size  from  less  than  ten  lines  of  code  to  many  thousands  of 
executable  statements.  Four  types  of  modules  exist  in  NASTRAN,  each  with  its 
specific  range  of  functions:  preface  modules;  executive  operation  modules;  utility 
modules;  and  analysis-oriented  modules  (Ref.  11).  Communication  between  modules 
is  performed  using  data  block  files  and  parameters  (stored  in  the  blank  common 
block).  Functional  modules  are  machine- independent,  and  comprise  about  90  percent 
of  the  code. 

Input/output  is  supported  through  GINO  (Generalized  Input/Output)  ,  a  compre¬ 
hensive  set  of  machine-dependent  utility  subroutines  developed  for  manipulating 
data  blocks  in  internal  NASTRAN  operations.  The  basic  unit  of  I/O  is  a  "logical 
record",  whose  word  length  is  variable.  Each  column  of  a  NASTRAN  matrix  data 
block  is  one  logical  record.  Two  subroutines  convert  special  NASTRAN  input  card 
formats  to  standard  FORTRAN  data  words  easily  handled  by  all  NASTRAN  input  pro¬ 
cessors.  All  NASTRAN  routines  must  use  GINO.  They  are  consequently  isolated 
from  the  actual  physical  hardware  and  such  concerns  as  blocking  factors  and 
device  characteristics.  Since  main  memory  is  used  as  a  scratch  pad  (no  module 
may  leave  values  in  main  memory) ,  GINO- formula ted  data  blocks  form  the  bulk  of 
the  intermodule  communication  (Ref.  14). 

NASTRAN  is  primarily  a  file-oriented  system,  using  mass  storage  devices  for 
nearly  all  major  data  transfers.  Four  types  of  input/output  systems  are  provided 
in  NASTRAN  (Ref.  11) :  external  user-supplied  interfaces;  internal  data  blocks 
(used  for  temporary  data  storage);  checkpoint/restart  files;  and  substructure 
operating  files  (storing  all  data  necessary  for  a  complete  multi-stage  sub¬ 
structuring  analysis).  NASTRAN  has  a  sophisticated  checkpoint  and  restart 
capability,  offering  the  user  four  general  types  of  restarts:  unmodified  restart, 
pseudo-mod ified  restart,  modified  restart  and  rigid  formal  switch.  The  reader 
is  referred  to  References  7,  8,  and  11  for  details  of  restarting. 

NASTRAN  is  not  a  core  program.  It  is  designed  to  dynamically  allocate 
core  so  that  at  execution  time,  all  the  core  memory  that  can  be  made  available 
will  be  used.  "Open  core"  memory  management  means:  a  contiguous  block  of 
randomly  addressable  working  storage  defined  by  a  labeled  common  block,  whose 
length  is  a  variable  determined  by  the  NASTRAN  executive  subroutine  KORSZ. 
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Main  memory  is  treated  as  a  large  single-dimensional  scratch  array  by  all 
modules.  The  length  is  communicated  dynamically  to  the  module  by  the 
Executive  System.  Fixed-dimension  statements  for  arrays  are  not  allowed 
in  any  module  as  all  user  input  is  open-ended  (Ref.  14).  Spill  logic  is 
provided  to  transfer  data  to  scratch  files  if  complete  core  allocation  is 
impossible.  NASTRAN  was  not  originally  designed  to  operate  under  a  time-sharing 
system. 

Each  module  in  NASTRAN  is  assigned  separate  overlay  segments  independent 
of  the  other  modules.  During  execution,  only  the  module  code  and  the  required 
utility  subroutines  are  loaded  into  core.  Furthermore,  a  common  block  assigned 
by  the  module  is  loaded  below  the  current  code  for  use  in  an  "open-core"  storage 
space.  Open  core  exists  from  this  point  to  the  maximum  core  space  requested 
for  the  execution.  All  matrix  and  table  data  required  for  a  given  module  oper¬ 
ation  are  brought  into  open  core  for  temporary  storage.  Results  are  moved  from 
open  core  to  data  block  files  for  access  by  other  modules.  A  typical  NASTRAN 
overlay  tree  is  illustrated  in  Figure  4-2.  Examples  of  modules  are  shown;  these 
are  divided  into  separate  parts,  consisting  of  a  small  "driver"  routine  which 
calls  the  major  sections  of  code  located  near  the  bottom  of  core.  This  scheme 
allows  the  large  and  complex  matrix  utility  routines  to  be  used  by  more  than  one 
module  in  a  link,  yet  still  allows  the  largest  possible  open  core  space  for  each 
module  (Ref.  11). 

COSMIC/NASTRAN  and  the  open  core  concept  are  implemented  on  three  major 
computer  systems;  IBM  360/370  series,  UNIVAC  1108/1110  series,  and  CDC  6600/ 
CYBER  series.  The  NASTRAN  program  is  divided  into  a  series  of  logical  pieces 
called  "links".  Each  link  is  a  complete  program  or  load  module,  and  contains 
its  own  overlay  structure  and  its  own  root  segment  (the  set  of  subprograms 
which  is  always  resident  in  main  memory  for  that  link).  Communication  between 
links  occurs  through  computer  files.  Table  4-1  illustrates  some  characteristics 
of  NASTRAN  implementation  on  the  three  computer  systems. 

Table  4-2  compares  some  program  architecture  features  of  NASTRAN  with 
four  general  purpose  codes:  MARC,  ADINA,  AGGIE,  and  ANSYS.  The  only  code  among 
these  give  to  offer  a  general  matrix  manipulation  capability  is  NASTRAN. 


Open  Core 
Part  2 


FIGURE  4-2.  ILLUSTRATION  OF  NAS TRAN  OVERLAY  SYSTEM. 
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Two  recent  articles  which  survey  and  compare  general  purpose  code  capabilities 
and  features  are  Fredriksson  and  Mackerle  (Ref.  21)  and  Dunder  and  Belonogoff 
(Ref.  22). 


TABLE  4-1 


NASTRAN  IMPLEMENTATION  ON  3  MAJOR  COMPUTER  SYSTEMS 


IBM  360/370 

Series 

UNIVAC  1108/1110 
Series  (Exec  8  OS) 

CDC  CYBER/6000 

Series  (NOS) 

Word  size 

32  bits 

36  bits 

60  bits  (single  precis 

Character 

capacity 

8  bits/character 

4  characters/ 
word 

6  bits/character 

6  characters/ 
word 

6  bits/character 

10  characters/word 

Number  of 
links 

16 

15 

15 

Overlay 

Open  core, 
labeled  common 
block 

Open  core 
segment  loader , 
common  block 

Segmentation  loader, 
segments  loaded 
dynamically 

Input/Output 

GINO  consists 
of  3  decks: 

.GINO 

.NASTIO 

.OPEN 

3  categories  of 

I/O  routines: 
.Identify  and  obtain 
logical  files 

GNFIAT 

.Perform  actual  I/O 
(GINO,  101108) 

. SGINO:  generate 
plot  tape 

4  ways: 

.Card  input,  printed/ 
punched  output, 
binary  output 
.All  other  I/O  uses 
GINO 

.SGINO:  Plot  output 
. SOFIO:  substructuring 
I/O 
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SECTION  5 

FUNCTIONAL  DESCRIPTION 

This  section  briefly  highlights  some  NASTRAN  functional  descriptions, 
in  order  to  give  the  reader  an  overview  of  the  versatility  and  capabilities 
of  the  code.  These  descriptions  include: 

•A  rigid  format  -  What  is  it?  How  many  are  there?  and  a  typical  DMAP 
sequence  for  linear  statics  analysis 
•Direct  Matrix  Abstraction  Program  -  What  is  DMAP?  and  a  short  example 
.Structural  Elements  -  a  summary  table  giving  element  categories,  names, 

descriptions,  number  of  nodes  and  degrees  of 
freedom  at  each  node,  material  property,  and  origin 
.Constitutive  library  -  material  property  options 
.Constraints  and  partitioning 
.Matrix  operations 

.Eigenvalue  extraction  methods,  compared  to  other  general 
.Time  integration  scheme  *  purpose  codes 

.Plotting  capabilities 
.Restart  capability 
.Nonlinear  capabilities 

.Special  NASTRAN  features  -  not  usually  found  in  other  codes 
.Substructuring  -  automated  multi-stage  substructuring  analysis 

5.1.  Rigid  Formats 

A  total  of  20  rigid  formats  exists  in  Level  17.5  of  COSMIC/NASTRAN.  A 

rigid  format  is  an  established  sequence  of  DMAP  instructions  stored  in  the 

Executive  System  to  perform  some  standardized  analyses,  such  as: 

.Static  analysis  -  conventional  linear  static  analysis;  with  inertia  relief 

(inertia  effects  of  unconstrained  rigid  body  accelerations) 
with  differential  stiffness  (nonlinear  effects  of  large 
deflections) 

.Normal  modes  analysis  (linear  elastic  models) 

.Buckling  analysis 

.Piecewise  linear  (material  plasticity)  analysis 

.Direct  complex  eigenvalue  analysis 

.Direct  frequency  and  random  response 

.Direct  transient  analysis 

.Modal  complex  eigenvalue  analysis 

.Modal  frequency  and  random  response 

.Modal  transient  analysis 

.Normal  modes  with  differential  stiffness 

.Static  analysis  using  cyclic  symmetry 

.Normal  modes  analysis  using  cyclic  symmetry 

.Heat  transfer  analysis  -  linear  or  nonlinear  steady-state;  transient 
.Aeroelasticity  analysis  -  modal  aerodynamic  flutter  analysis;  modal 

aeroelastic  response  due  to  applied  loads  and  gusts 
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Rigid  formats  are  identified  with  an  integer  number  and  an  analysis  approach. 
There  are  three  approaches:  15  rigid  formats  in  displacement  (APProach  DISP) ; 

3  in  heat  transfer  (APProach  HEAT) j  and  2  in  aerodynamic  (Approach  AERO) .  Page 
4.12-3  of  Reference  11  contains  an  excellent  table  which  summarizes  the  various 
modeling  options  offered  versus  rigid  format  number.  These  rigid  formats  provide 
the  typical  user  with  a  tremendous  range  of  analysis  options,  using  pre-established 
module  call  sequences  (the  usual  case  for  all  but  the  expert  user).  If  he  chooses, 
the  expert/experienced  user  can  modify  the  rigid  formats  by  DMAP  alters  for  a 
particular  problem. 

Table  5-1  shows  the  DMAP  sequence  for  a  typical  rigid  format,  in  this  case 
the  flow  sequence  of  the  30  major  functional  modules  for  linear  statics  analysis. 
Note  that  the  modules  have  been  classified  into  pre-processing,  analysis,  and 
post-processing  modules.  This  is  the  only  rigid  format  with  a  fully-stressed 
design  optimization  capability,  and  the  optimization  loop  is  shown.  SSG3  is  the 
key  analysis  module  here  since  it  solves  for  the  independent  displacements.  The 
rigid  format  for  normal  modes  analysis  has  a  very  similar  DMAP  sequence  to  that 
of  linear  statics.  The  key  analysis  module  for  normal  modes  analysis  is  named 
READ  (real  eigenvalue  analysis  -  displacement),  consisting  of  some  52  subroutines. 
5.2.  Direct  Matrix  Abstraction  Program  (DMAP) 

DMAP  is  one  of  the  most  powerful  matrix  manipulation  tools  offered  in  any 
code.  It  is  the  programming  language  of  NASTRAN.  DMAP  is  simply  a  language  of 
macro  instructions  which  enables  the  user  to  sequence  any  combination  of  matrix 
operations.  A  DMAP  instruction  has  the  form: 


(Module  Name)  (Input  Data  Blocks) /(Output  Data  Blocks) /(Parameter  List)? 

where  the  slash  is  used  as  a  delimiter  between  data  blocks  and  the  dollar  sign 
is  used  to  terminate  the  DMAP  instruction. 

A  short  DMAP  example  is  explained  below  (Ref.  14): 

To  compute  [c]  =  [X]  +  [B] 

[D]  -  [A]  [C] 

The  DMAP  sequence  is: 

BEGIN  $  START  DMAP 

ADD  A,B/C/$  ADDS  A  TO  B 


OPTIMIZATION  LOOP 
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TABLE  5-1.  DMAP  SEQUENCE  FOR  RIGID  FORMAT  1  -  LINEAR  STATICS 


Module 

Name  _  _ FUNCTION 


GP1 


Pre¬ 

processing 

Modules 


3P2 
?LOT 
IP  3 
TA1 

2£TPR1 


Coordinate  system,  grid  point  locations,  relate  internal 
to  external  grid  points 
Element  connection  table 
Undeformed  structure  plots 
Static  loads,  grid  point  temperatures 

Element  tables  for  use  in  matrix  assembly  and  stress  recovery 
Phase  1  property  optimization,  initialization  check 


SMG 

SMA 


Analysis 

Modules 


DMAP 

Loop 


3PWG 

5MA3 


IP  4 

3PSP 

4CE1 

4CE2 

SCE1 

3MP1 

RBMG1 

XBMG2 

RBMG3 

3SG1 

3SG2 


3SG3 

3DR1 

3DFDR 

3DR2 


Element  stiffness  and  mass  matrices,  for  later  assembly 

Element  matrix  assembler  assembles  stiffness  matrix,  grid- 

point  table  singularity 

Weight  and  balance  information 

Adds  general  elements  to  stiffness  matrix 

Displacement  sets,  MPC  equations,  enforced  displacement  vector 
Determines  if  possible  grid  point  singularities  remain 
Positions  MPC  equations 
Partitions  stiffness  matrix 
Partitions  out  SPCs 

Partitions  constrained  stiffness  matrix 
Partitions  out  free-body  supports 
Decomposes  constrained  stiffness  matrix 

Rigid  body  transformation  and  check  matrices,  error  ratio 
Generates  static  load  vectors 

Applies  constraints  to  static  load  vectors,  calculates 
determinate  reactions 
Solves  for  independent  displacements 
Recovers  dependent  displacements  and  SPC  forces 

Calculates  grid  point  force  balance  element  strain  energy  as 
requested 

Calculates  element  forces  and  stresses 


§DR3 


Post¬ 

processing 

Modules 


XYTRAN 

XYPLOT 

DPTPR2 

PLOT 


Prepares  requested  output  sorted  by  grid  point  number  of 
element  number 

Prepares  input  for  X-Y  plots 

Prepares  X-Y  plots  of  displacements,  stresses,  forces,  SPC 
forces  vs.  subcase 

Performs  Phase  2  property  optimization 
Deformed  structure  and  contour  plots 


J 


a* 
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MATPRN  C//$  PRINT  C 

MPYAD  A/C ,/D/$  COMPUTE  D 

MATPRN  D//$  PRINT  D 

END  $  TERMINATE 

This  simple  DMAP  example  provides  an  illustration  of  the  use  of  DMAP  which  is 
easy  to  understand  even  for  the  relatively  inexperienced  novice  and  common  users 
of  NASTRAN.  In  addition  to  matrix  arithmetic,  DMAP  in  NASTRAN  can  be  used  for 
tasks  categorized  as:  executive,  utility,  and  structural  (e.g.,  assemble  tables, 
calculate  element  data) . 

A  rigid  format  sequence  may  be  modified  by  the  user  using  the  ALTER  capability. 
A  library  of  useful  pre-defined  DMAP  alter  sequences  exists,  and  these  sequences 
are  called  RFALTERS.  The  user  must  merge  RFALTERS  into  the  Executive  Control  Deck 
by  using  a  machine-dependent  operating  system  utility  for  merging  files. 

5.3.  Structural  Elements 

Table  5.2  gives  a  summary  of  the  structural  elements  available  in  Level  17.5 
of  COSMIC/NASTRAN.  Of  the  33  elements  listed,  4  are  higher-order  elements  (TRIM6, 
TRPLT1 ,  TRSHL,  TORDRG)  and  4  are  isoparametric  elements  (QDMEM1 ,  IHEX1,  IHEX2 , 

IHEX3) .  Nineteen  elements  have  heat  transfer  capability.  For  simplicity,  fluid 
elements  and  some  scalar  elements  have  not  been  included  in  Table  5-2,  such  as: 
viscous  damper,  scalar  spring  and  mass,  concentrated  mass,  general  element,  and 
dummy  element.  By  today's  standards,  the  element  library  in  NASTRAN  is  rated  only 
fair.  Many  of  the  elements  are  about  10  to  15  years  old  and  have  existed  in 
NASTRAN  since  its  inception.  One  glaring  weakness  is  the  absence  of  a  quadrilateral 
solid-of-revolution  element  for  plane  strain,  plane  stress,  or  axisymmetric 
analyses.  The  trapezoidal  ring  element  TRAPRG,  which  requires  two  sides  of  the 
element  to  be  parallel,  is  simply  too  restrictive  and  not  versatile  enough.  A 
major  portion  of  this  study  was  to  perform  advanced  evaluation  exercises  on 
selected  NASTRAN  elements.  The  elements  selected  for  advanced  evaluation  are 
noted  in  Table  5-2,  and  the  reader  is  referred  to  the  appropriate  subsection  in 
Section  7  for  evaluation  results.  Page  3.2-2  of  Reference  11  contains  a  detailed 
summary  table  of  NASTRAN  finite  element  characteristics  (e.g.  thermal  loads,  lumped 
or  consistent  mass,  differential  stiffness,  plasticity,  heat  transfer,  etc.).  A 
library  of  more  than  60  different  finite  element  formulations  is  offered  in 
COSMIC/NASTRAN  Level  17.5. 
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TABLE  5-2.  STRUCTURAL  ELEMENTS  SUMMARY 
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5.4.  Constitutive  Library 

As  already  indicated  in  Table  5-2,  NASTRAN  allows  three  basic  types  of 
material  properties:  isotropic;  orthotropic  (for  TORDRG  and  solid-of-revolution 
elements  only);  and  anisotropic  (2-D  flat  surface  elements  only).  In  addition 
to  these  elastic  properties,  NASTRAN  material  data  cards  are  also  used  to  define 
mass  density,  thermal  expansion  coefficients,  and  stress  limits.  Each  material 
is  given  an  unique  identification  number  which  may  be  referenced  by  any  number 
of  property  cards.  Material  data  cards  are  identified  in  the  Bulk  Data  Deck 
as  follows: 

.Temperature-independent  material  properties  (MATi) 

MAT1  elastic  isotropic 

MAT2  anisotropic,  2-D  flat  surface  elements  only 

MAT3  elastic  orthotropic,  TORDRG  and  solid-of-revolution  elements  only 
MAT4  heat  transfer  analysis,  isotropic  properties 
MAT5  heat  transfer  analysis,  anisotropic  properties 

.Temperature-dependent  material  properties 

MATTi  temperature  dependence  in  conjunction  with  MATi 
TABLEMi  tabular  functions  of  temperature-dependent  properties 

.Stress-dependent  material  properties 

MATSi  stress  dependence  in  conjunction  with  MATi 

5.5.  Constraints  and  Partitioning 

NASTRAN  offers  the  user  a  variety  of  options  to  constrain  and  reduce  his 
problem.  These  constraints  are  used  to  specify  fixed  boundary  conditions,  to 
eliminate  matrix  singularities,  to  define  rigid  elements,  and  to  support  free 
body  motion.  The  most  important  of  these  are:  multiple  point  constraints  (MPC's), 
single  point  constraints  (SPC's),  OMIT  and  ASET.  A  multipoint  constraint  is  a 
linear  relationship  between  two  or  more  displacement  degrees  of  freedom.  Closely 
related  to  MPC's  are  the  rigid  elements  (RIGDR,  RIGD1-3)  in  NASTRAN,  which  are 
very  stiff  connections.  These  can  be  used  to  model  levers,  pulleys,  gear  trains, 
and  rigid  links  to  remove  matrix  ill-conditioning.  SPC's  are  vectors  of  enforced 
displacements,  any  or  all  of  whose  elements  may  be  zero.  SPC's  commonly  represent 
structural  boundary  conditions  of  zero  displacement  and  slope,  or  enforced  defor¬ 
mations  at  nodes. 

For  matrix  reduction  the  user  can  select  either  the  OMIT  card  or  the  ASET 
card,  but  not  both.  OMIT  specifies  degrees  of  freedom  to  be  reduced  out  of 


5-8 


the  analysis  set,  while  ASET  specified  independent  degrees  of  freedom  to  be 
retained.  Since  they  represent  complementary  sets,  the  user  should  specify 
the  smaller  of  the  two.  The  specified  degrees  of  freedom  for  elimination  by 
reduction  are  restricted  to  only  nonconstrained  coordinates.  Another  bulk  data 
card,  SUPORT,  may  also  be  used  to  specify  the  degrees  of  freedom  which  will 
remove  rigid  body  motion.  SUPORT  is  used  to  supply  temporary,  nonredundant 
constraints  on  free  body  motion.  It  is  required  for  inertial  relief  analysis 
(Rigid  Format  2) ,  and  recommended  for  use  in  modal  analysis  of  free  bodies. 

Errors  in  constraint  data  are  not  automatically  corrected  by  NASTRAN.  Con¬ 
flicting  constraint  data,  when  detected,  are  treated  as  fatal  errors  (Ref.  11). 
Grid  point  singularity  tests  are  performed  to  detect  null  columns  in  the  stiffness 
matrix  resulting  from  deficient  constraints  or  element  connectivities. 

5.6.  Matrix  Operations 

An  unique  feature  of  NASTRAN  is  the  user's  option  to  manipulate  matrices 
using  DMAP.  Therefore,  theoretically,  he  can  solve  any  numerical  analysis  problem 
which  can  be  expressed  in  matrix  form.  Elements  of  matrices  can  be  specified  by 
the  DMI  card.  The  user  can  build  a  DMAP  sequence  using  any  combination  of  43 
matrix  handling  subroutines.  Design  philosophy  of  the  modules  is  strictly  based 
on  handling  large  matrices.  All  matrices  are  stored  on  peripheral  devices  (tapes, 
disks  and/or  drums)  by  columns,  and  are  packed  in  non-zero  strings.  Matrix 
sparsity  and  bandedness  are  emphasized  and  utilized.  The  user  may  select  single 
or  double  precision  control  by  the  executive  control  card  PREC;  this  is  the  only 
method  offered  in  NASTRAN  to  combat  round-off  error  accumulation.  This  matrix 
handling  feature  in  NASTRAN  is  elaborate,  sophisticated,  and  powerful.  (The  only 
other  general  purpose  code  which  offers  a  similar  matrix  operations  capability 
is  SPAR.) 

5.7.  Eigenvalue  Extraction  and  Time  Integration 

Table  5-3  compares  the  eigenvalue  extraction  methods  and  direct  time  inte¬ 
gration  scheme  of  NASTRAN  with  eight  other  U.  S.  general  purpose  codes.  This 
comparison  reveals  that  while  NASTRAN  offers  more  eigenvalue  extraction  routines 
(4)  than  any  other  code,  its  one  time  integration  scheme  is  less  than  an  average 
of  two  schemes  offered  by  the  other  codes.  The  most  recent  method  of  eigenvalue 
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extraction,  subspace  iteration  (Ref.  51),  is  implemented  in  only  three  codes: 

ADINA,  SAP6 ,  and  EASE2. 

For  direct  time  integration,  three  methods  have  gradually  emerged  to  be  the 
most  popular  among  U.  S.  code  developers:  Newmark  Beta,  Wilson  and  Houbolt. 

Newmark  Beta  and  Wilson  are  "implicit"  methods,  where  a  matrix  system  is  solved, 
one  or  more  times  per  step,  to  advance  the  solution.  Both  algorithms  are  un¬ 
conditionally  stable,  but  recent  research  has  indicated  that  the  Wilson  algorithm 
introduces  damping  into  the  solution,  is  sensitive  to  time  step  size,  and  may 
require  more  steps  to  obtain  the  same  accuracy  as  the  Newmark  method.  The  Houbolt 
method  is  a  third-order  backward  difference  scheme  which  is  unconditionally  stable 
for  all  time  step  sizes.  However,  it  introduces  damping  into  the  solution,  in 
addition  to  having  a  tendency  to  remove  higher  modes  from  the  system.  The  central 
difference  scheme  offered  in  ADINA  and  MARC  is  an  "explicit"  method,  where  the 
solution  may  be  advanced  without  storing  a  matrix  or  solving  a  system  of  equations. 
It  requires  very  small  time  steps,  but  can  be  cheaper  in  cost  per  time  step.  For 
nonlinear  dynamic  problems,  the  most  sophisticated  codes  are  ADINA,  MARC,  and  a 
new  code  called  ABACUS  being  developed  by  Hibbitt  and  Karlsson,  Inc. 

Table  5-4  is  a  detailed  comparison  of  the  four  NASTRAN  eigenvalue  extraction 
methods.  The  method  desired  for  a  particular  solution  is  selected  by  the  user  on 
the  EIGR  card.  The  two  tridiagonalization  methods  (Givens  and  FEER)  are  cate¬ 
gorized  as  "transformation"  methods  (obtaining  all  the  eigenvalues  at  once)  and 
can  be  used  only  for  real  matrices.  The  newer  FEER  method  is  a  tridiagonalization 
procedure  based  on  the  Rayleigh-Ritz  method.  It  combines  the  best  features  of  the 
inverse  power  and  Givens  methods.  The  other  two  methods  (inverse  power  and  deter¬ 
minant)  are  "tracking"  schemes  and  can  also  be  applied  to  complex  matrices.  Proper 
selection  of  the  extraction  method  depends  on  the  number  of  eigenvalues  desired, 
matrix  character,  problem  size,  and  accuracy  demanded;  Table  5-4  offers  some 
hints  in  this  choice.  This  subject  is  discussed  in  more  detail  in  Ref.  51  and 
Chapter  13  of  Ref.  13. 
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TABLE  5-3.  COMPARISON  OF  EIGENVALUE  EXTRACTION  AND  TIME  INTEGRATION  METHODS 
IN  TEN  U.  S.  GENERAL  PURPOSE  CODES 


CODE 

EIGENVALUE  EXTRACTION  METHODS 

TIME  INTEGRATION  SCHEMES 

1 

NASTRAN 

1. 

Tridiagonal  (Givens) 

Newmark  Beta  (6  =  7)  method 

(Ref.  7,11) 

2. 

Inverse  Power  with  Shifts 

3. 

Determinant  Method 

.For  heat  transfer  transient 

4. 

Tridiagonal  Reduction 

analysis:  user’s  choice 

Method  -  or  Fast  Eigen- 

0<B<1 

value  Extraction  Routine 

(FEER) 

2 

MARC 

Inverse  power  sweep 

1. 

Central  difference  . 

(Ref.  15) 

2. 

Newmark  Beta  (g  =  — ) 

3. 

Houbolt  method 

3 

ADINA 

1. 

Determinant  Search  Method 

1. 

Central  difference 

(Ref.  17,35) 

(with  Sturm  sequence) 

2. 

Wilson  0=1.4 

2. 

Subspace  iteration 

3. 

Newmark  Beta 

4 

AGGIE  I 

(unnamed)  "tracking"  scheme 

1. 

Newmark  Beta 

(Ref.  19) 

2. 

Wilson  0  *  1.4 

5 

ANSYS 

Jacobi 

Houbolt 

(Ref.  20) 


6 

STARDYNE 

1. 

Householder-QR  Method 

None  (planned  for  implemen- 

(Ref.  23) 

2. 

(tr idiagonalization) 

Inverse  iteration  (same  as 
NASTRAN' s  inverse  power  with 
shifts) 

tat ion,  but  undocumented) 

3. 

Lanczos  modal  extraction  method 
(same  as  NASTRAN’ s  FEER  method) 

7 

SAP  6 

1. 

Determinant  Search 

Wilson  0  =  1.4 

(Ref.  24) 

2. 

Subspace  Iteration 

8 

mm 

H 

Determinant  Search  l. 

Newmark  Beta 

2. 

Subspace  Iteration  2. 

Wilson  0 

9 

SDRC/SUPERB 

1. 

Determinant  Tracking  (with 

None  (NASTRAN  &  ANSYS  inter- 

(Ref.  26) 

2. 

Sturm  sequence  property) 

Jacobi  (with  Guyan  reduction) 

face  packages  offered) 

10 

MSC/NASTRAN 

1. 

Modified  Givens  (MGIV) 

Newmark  Beta  Method 

(Ref.  41) 

2. 

Inverse 

Analytical  solution  of 

3. 

Givens 

equation  (Quasi-closed  form) 
(also  has  component  mode  syn 
thesis  and  superelement  dy¬ 
namic  capability  and  general 
ized  dynamic  reduction 

TABLE  5-4.  COMPARISON  OF  NASTRAN  EIGENVALUE  EXTRACTION  METHODS 
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5.8.  Plotting  Capabilities 

The  user  specified  undeformed  or  deformed  structure  plots  at  the  end  of 
the  Case  Control  Deck.  The  structure  plot  request  packet  begins  with  the  card 
OUTPUT  (PLOT).  The  plotting  capabilities  in  NASTRAN  are  excellent,  comprehensive, 
and  well-documented  (Ref.  7,  8,  11)  with  many  sample  plot  card  setups.  The 
following  types  of  plots  may  be  selected: 

.Undeformed  geometry  plots  for  all  rigid  formats 

.Static  deformations,  and  stress/displacement  contours 

.Modal  deformations  (eigenvectors  or  mode  shapes) 

.Transient  response  or  frequency  response-vectors  or  deformed  shape 
for  specified  times  or  frequencies 

.X-Y  graphs  of  transient/frequency  response 

.V-f  and  V-g  graphs  for  flutter  analysis 

.Topological  displays  of  matr ices-showing  nonzero  element  locations 
when  the  SEEMAT  utility  module  is  requested 

NASTRAN  plotting  is  compatible  with  most  plotter  hardware  like  Stromberg- 

Carlson  and  CALCOMP.  The  user  has  a  choice  of  three  types  of  projections: 

orthographic,  perspective  and  stereoscopic  (microfilm  plotters  only).  Plot 

labeling  and  scaling  are  user-controlled  or  automatically  provided.  NASTRAN 

plotting  does  not  have  hidden-line  capability,  which  exists  in  codes  such  as 

MARC  and  ANSYS.  Currently,  no  interactive  pre-processing  and  post-processing 

capability  exist  (see  Sec.  7.9)  or  are  planned,  and  NASTRAN  runs  as  well  as 

plots  are  used  in  the  batch  mode.  Cost  estimates  for  plotting  are  missing  from 

the  documentation;  plotting  costs  can  be  considerable. 

5.9.  Restart  Capability 

NASTRAN  contains  a  sophisticated  restart  capability.  Restarts  are  effective 
for:  (1)  continuing  problem  execution  having  an  unscheduled  interruption  caused 
by  a  data  error;  (2)  requesting  additional  information  for  a  problem  already  com¬ 
pleted;  (3)  running  additional  load  cases  in  static  analysis;  and  (4)  extracting 
real  eigenvalues  for  additional  frequency  ranges  in  normal  modes  analysis.  Re¬ 
starts  are  ineffective  when  the  problem  is  small,  or  when  changes  are  made  in  the 
element  properties  or  grid  point  information.  A  CHKPNT  YES  card  in  the  Executive 
Control  Deck  will  save  information  on  a  new  problem  tape  NPTP,  subsequently  renamed 
an  old  problem  tape  OPTP  in  a  future  run.  One  disadvantage  of  using  restarts 
is  the  necessity  of  maintaining  checkpoint  files. 


5.10.  Nonlinear  Capabilities 


The  nonlinear  capabilities  in  NASTRAN  are  limited.  A  rough  guess  is  that 
probably  less  than  5  percent  of  all  NASTRAN  users  has  ever  attempted  the  non¬ 
linear  usage.  Two  types  of  static  nonlinearities  are  offered:  differential 
stiffness  and  material  plasticity.  The  geometric  stiffness  effects  treated  by 
Rigid  Formats  4  and  5  provide  the  user  with  a  second  order  approximation  to 
the  nonlinear  effects  of  large  deflection.  Rigid  Format  6  offers  a  piecewise 
linear  analysis  capability  (material  plasticity).  Only  isotropic  materials  are 
allowed,  and  the  capability  is  restricted  to  certain  elements:  ROD,  TUBE,  BAR, 
and  plate  elements  (based  on  2-D  plasticity  theory  by  J.  L.  Swedlow) .  The 
material  properties  are  assumed  to  be  stress  dependent.  The  stiffness  matrix 
is  assumed  constant  over  each  load  increment.  After  each  increment,  a  new  stiff¬ 
ness  matrix  is  generated  based  on  the  current  state  of  stress  in  each  element. 

All  static  load  options  are  allowed,  except  temperature  and  enforced  element 
deformations. 

For  transient  analysis,  the  module  TRD  provides  for  four  types  of  nonlinear 
elements,  mainly  for  control  systems  simulation.  Nonlinear  effects  are  treated 
as  an  additional  applied  load  vector,  whose  elements  are  functions  of  either 
displacements  or  velocities.  One  additional  nonlinear  capability  is  nonlinear 
steady-state  heat  transfer  analysis,  which  allows  for:  temperature-dependent 
conductivities  for  the  elements,  nonlinear  radiation  exchange,  and  a  limited 
use  of  multipoint  constraints.  Nonlinear  elements  such  as  gaps  and  friction 
elements  do  not  exist  in  COSMIC/NASTRAN .  The  NASTRAN  nonlinear  capabilties 
offered  are  primitive  when  compared  to  those  of  current  nonlinear  codes  such  as 
ADINA,  AGGIE,  ABACUS,  ANSYS,  and  MARC. 

5.11.  Special  NASTRAN  Features 

NASTRAN  contains  certain  special  or  unique  features  not  commonly  found  in 
most  general  purpose  codes.  These  are  briefly  surveyed  in  Table  5-5.  In  general, 
the  most  likely  users  of  these  features  are  the  expert/experienced  ones.  The 
common  and  novice  users  will  probably  find  most  of  these  features  too  difficult  to 
use.  These  features  represent  a  considerable  investment  in  time  and  money,  and 
are  no  doubt  valuable  to  certain  small  pockets  of  the  user  community.  However, 
their  overall  effectiveness  and  extent  of  usage  are  difficult  to  assess,  and  were 
not  evaluated  in  the  study. 


5.12.  Substructurim 


The  automated  multi-stage  substructuring  system  (Ref.  7,  11,  27)  in  NASTRAN, 
developed  by  Universal  Analytics,  Inc.,  is  very  powerful  yet  flexible  and  user- 
oriented.  It  incorporates  an  automated  multi-stage  modal  synthesis  procedure 
(Ref.  28)  for  the  dynamic  analysis  of  very  large  models,  and  a  multi-stage 
component  mode  synthesis  method  (Ref.  29).  Table  5-6  compares  the  NASTRAN  auto¬ 
mated  multi-stage  substructuring  capability  versus  that  offered  by  ANSYS.  Both 
codes  have  highly-developed  substructuring  capabilities  which  have  been  used 
widely.  Similarities  and  differences  are  seen  in  the  design  philosophy  of  the 
two  codes  and  both  advantages  and  disadvantages  are  noted.  Substructuring  is 
suggested  only  for  the  expert/experienced  users.  A  learning  curve  is  invariably 
necessary.  However,  for  certain  classes  of  problems,  substructuring  offers 
tremendous  benefits  and  savings. 
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TABLE  5-5.  SPECIAL  NAS TRAN  FEATURES 


FEATURE 

DESCRIPTION 

1. 

Static  analysis 
with  inertia  relief 

.Treats  structural  models  not  fully  constrained. 

.Computes  inertia  effects  of  unconstrained  rigid  body 
accelerations 

2. 

Manual  single-stage 
substructure  analysis 

.Phase  I  analyze  each  substructure,  establish  master 

degrees  of  freedom 

.Phase  II  analyze  "pseudo-structure" 

.Phase  III  back  substitute  for  substructure  displace¬ 
ments  and  stresses 

•Can  be  used  for  all  rigid  formats  except  piecewise 
linear  analysis 

3. 

Fully  stressed  design 
optimization  (Rigid 
Format  1  -  Linear 
Statics) 

(Ref.  11  Sec.  26) 

.Very  simple  design  algorithm  to  resize  elements 
.For  each  design  iteration,  change  each  element  cross- 
sectional  properties  to  get  a  limit  stress  somewhere 
within  element 

.Procedure  works  reasonably  well  for  statically  deter¬ 
minate  structures,  acceptable  for  slightly  indeter¬ 
minate,  and  poor  for  indeterminate  structure  with  high 
redundancy 

4. 

Cyclic  Symmetry 

(Ref.  7,11,40) 

.Linear  problems  only;  static  and  normal  modes  analyses 
.Uses  finite  Fourier  transformation 

.Symmetric  with  respect  to  axis:  rotational,  dihedral 
symmetry 

.Model  1  segment;  use  uncoupled,  transformed  equations 

5. 

Component  Mode 
Synthesis 

(Ref.  29) 

.Design  tool  -  evaluate  critical  component  modes  in 
total  response 

•Divide  structure  into  separate  components,  reduce  order 
of  component  matrices,  combine  them  using  compatibility 
at  common  boundary  points 

6. 

Automated  multi-stage 
substructure  analysis 

(Ref.  11,27,28) 

.3  phases  like  "2",  substructure  operating  file  (SOF) 
created 

.User  can  repeatedly  combine  and  reduce  structures 
.Linear  static  and  dynamic  analyses 

.Sequence  control  now  automated  using  Substructure  Case 
Control  deck 

7. 

Representation  of 
part  of  a  structure 
by  its  vibration 
modes 

.Use  modal  information  from  other  analyses  or  vibration 
test  data 

.3  cases:  all  connection  coordinates  free,  restrained, 
or  mixed 

.Restrictions:  linearity,  conservation  of  energy, 
reciprocity 
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TABLE  5-5  cont'd. 


FEATURE 

DESCRIPTION 

8.  Representation  of 
control  systems 

.Linear  and  nonlinear  control  systems  simulated 
•Properties  of  control  system  treated  in  quadratic  format 
like  in  dynamics:  Mp*  +  Bp  +  K 

9.  Structure/fluid 
interaction 
(Hydroelastic 
analysis) 

.Compressible  fluids  in  axisymmetric  tanks 
.Compressible  fluids  in  rotationally  symmetrical  cavities 
with  slots 

10.  Acoustic  cavity 
analysis 

.For  each  element,  pressure  field  assumed  to  vary  linearly 
over  the  cross  section  and  sinusoidally  around  the  axis 
in  circumferential  direction. 

.Slot  portion  of  cavities  limited  to  certain  shapes 

11.  Aeroelastic 
analysis 

•Can  be  used  in  conjunction  with  structural  analysis 
.2  rigid  formats:  modal  flutter  analysis  (by  3  methods); 
modal  transient  and  frequency  response  analysis  of 
aeroelastic  models 

12.  Special  elements 

.General  element  GENEL  .Aeroelastic  elements  AEROi 

.Dummy  elements  CDUMi  .Acoustic  elements  AXIFi 

•Rigid  elements  RIGDi  .Hydroelastic  elements  FLUID i 
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SECTION  6 

VERIFICATION  EXERCISES 

In  the  course  of  reviewing  the  copious  NASTRAN  documentation,  the  authors 
were  very  impressed  with  the  complexity  of  the  problems  solved  using  NASTRAN. 
However,  it  is  often  difficult  to  determine  the  accuracy  of  such  complex  solu¬ 
tions  as  there  are  often  no  theoretical  solutions  available.  Experimental 
results  are  limited  and  it  is  difficult  to  obtain  the  idealized  loading  and 
boundary  conditions  assumed  in  analysis.  Therefore,  a  limited  number  of 
structural  verification  problems  were  run.  The  problems  selected  could  be 

1 k 

easily  verified  against  "textbook"  solutions.  These  examples  are  included 
in  this  section. 

A  number  of  interesting  observations  can  be  made  based  upon  these  results. 
Example  1  indicates  that  the  accuracy  of  the  code  is  in  good  agreement  with 
the  closed  form  solution  for  this  problem.  Example  2  illustrates  the  effect 
that  MPC's  can  have  on  the  solution  for  certain  elements.  Example  3  was 
chosen  to  exercise  the  various  options  available  in  NASTRAN  for  eigenvalue 
extraction.  All  were  in  good  agreement  when  a  consistent  mass  matrix  was  used. 


* 

The  solutions  to  these  problems  using  the  ANSYS  general  purpose  computer  program 
are  available  from  the  open  literature  and  are  included  for  comparison. 


Type: 

Reference : 


EXAMPLE  1.  BENDING  OF  A  TAPERED  PLATE 


Problem: 


Given: 


Static  analysis,  plate  elements 

Harris,  C.O. ,  Introduction  to  Stress  Analysis,  The  MacMillan  Co., 

New  York,  1959,  Page  114,  Problem  61. 

A  tapered  cantilever  plate  of  rectangular  cross-section  is  subjected 
to  a  load  P  at  its  tip.  Find  the  maximum  deflection  5  and  the  stress 
0x  in  the  plate. 

E  =  30  x  10^  psi,  L  =  20  in.,  b  =  3  in.,  d  =  0.5  in.,  P  =  10  lb. 


Results: 


<5  (in.) 

ax  (psi) 

Theory 

-0.0426666 

1600. 

ANSYS  (STIF  6) 

-0.0426668 

1600. 

(Ref.  30,  Problem  34) 

NASTRAN 

-0.0426672 

1625. 

(  TRIA2) 
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EXAMPLE  2.  ELONGATION  OF  A  SOLID  BAR 


Type:  Static  analysis,  solid  elements 

Reference;  Harris,  C.O. ,  Introduction  to  Stress  Analysis,  The  MacMillan  Co., 
New  York,  1959,  Page  237,  Problem  4. 


Problem: 


Given: 


A  tapered  bar  of  square  cross-section,  made  of  an  aluminum  alloy  for 
which  E  is  10,  400,000  psi,  is  suspended  from  a  ceiling.  An  axial 
load  Q  is  applied  to  the  end  of  the  bar.  Determine  the  maximum  axial 
deflection  6  in  the  bar  and  the  axial  stress  a at  mid-length. 

L  =  10  in.,  d  =  2  in.,  Q  =  10,000  lb. 


Problem  Sketch 

Finite 

Element  Model 

<5  (in.) 

ay  (psi) 

No.  of  input  cards 

Theory 

0.00480769 

4444. 

- 

ANSYS  (Ref.  30,  Problem 
37)  (STIF  5)  Difference 

0.00478878 

0.39% 

4441. 

0.067% 

33 
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EXAMPLE  2 
Results  cont'd. 

6 (in.)  Oy(psi)  No.  of  input  cards 


NASTRAN:  HEXA1  with  MPC's 
Difference 

Without  MPC's  to  enforce 
symmetry: 

1.  HEXAl  (Nodes  28,68) 

(Nodes  6,48) 
Difference 


0.00479047 

0.36% 


0.006179 
0.007493 
30 . 0%/56. 4% 


4441.  91 

0.067% 


4441.  67 


.067% 


2 .  1HEX1 


0. 00479047 


4441 


67 
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EXAMPLE  3.  NATURAL  FREQUENCIES  OF  A  CANTILEVER  BEAM 

Type :  Normal  mode  analysis  (eigenvalue  extraction) ,  beam  elements 

Reference:  Thomson,  W.  T.  Vibration  Theory  and  Applications,  Prentice-Hall,  Inc. 
Englewood  Cliffs,  N.  J. ,  1965,  Pages  275  and  357. 

Problem:  Determine  the  first  three  natural  frequencies  f^  of  a  uniform  beam 
clamped  at  one  end  and  free  at  the  other  end. 

Given:  E  =  30  x  106  psi,  I  =  1.3333  in4,  A  =  4  in2,  h  =  2  in. ,  L  =  80  in. , 
w  *  1.124  lb/in. 


i 


PROBLEM  SKETCH 


A 


4 

I 


(D  t  ©  f  ©  ♦  ©  *  ©  t  ©  I 


z  3  A  5 

Finite  Element  Model 


Dynamic  Degree 
of  Freedom 


(Uses  6  lateral  dynamic  degrees  of  freedom) 


Results: 


f  x  (Hz) 

f2(Hz> 

f 3 (Hz ) 

Theory 

10.247 

64.221 

179.82 

ANSYS  (Ref.  30,  Problem  51) 

10.247 

64.197 

180.14 

(STIF  3)  Difference 

None 

0.037% 

0.178% 

NASTRAN  (BAR) 

1.  Lumped  mass  model 

10.105 

61.437 

167.542 

2.  Inverse  power  method- 
consistent  mass 

10.203 

64.145 

179.905 

3.  Determinant  method 

10.233 

64.148 

179.905 

4.  Givens  method 

10.233 

64.148 

179.905 

5.  Unsymmetr ical  inverse  power 

10.233 

64.148 

179.905 

6.  Unsymmetr ical  determinant 

10.233 

64.148 

179.905 

SECTION  7 


ADVANCED  EVALUATION 

The  advanced  evaluation  consisted  of  a  comprehensive  literature  survey 
on  a  selected  group  of  NASTRAN  finite  elements  and  performing  a  series  of 
computer  runs  aimed  at  assessing  element  quality,  convergence,  completeness, 
and  efficiency.  As  much  as  possible,  the  framework  of  this  evaluation  is 
based  on  comparisons  with  other  well-known  benchmark  tests,  and  on 
recent  publications  on  simple  element  evaluation  tests  and  criteria 
(References  31  to  35) .  The  intent  of  this  section  is  to  document  in  one 
place  what  is  currently  known  about  the  quality  of  some  representative  COSMIC/ 
NASTRAN  structural  elements: 

QDPLT ,  TRSHL ,  TRIM6 ,  TRPLTl,  CONEAX,  TRAPRG,  TRAPAX,  HEXA1,  HEXA2 ,  IHEXl  and 

IHEX2. 

After  these  element  evaluation  results  are  described,  we  will  also  briefly 
discuss  current  minicomputer  versions  of  NASTRAN,  commercially  available  pre- 
and  post-processors  for  NASTRAN,  and  a  study  that  was  conducted  on  trends  of 
run  cost  and  CPU  time  versus  mesh  sizes. 

7.1.  Advanced  Evaluation  Tests 

To  aid  in  the  overall  assessment  of  any  finite  element,  Robinson  has 
suggested  several  evaluation  tests  and  the  following  assessment  points  (Ref.  31) 
before  "the  element  is  turned  loose  in  the  commercial  arena": 

1.  Element  shape. 

2.  Basic  assumption  (stress,  displacement,  or  mixed) . 

3.  Number  of  nodes. 

4.  Degrees  of  freedom,  types  of  freedom,  and  number  of  freedoms  per  node. 

5.  Are  the  freedoms  conventional  and  understandable? 

6.  Does  the  element  pass  the  standard  and  recognized  tests? 

7.  How  general  is  the  element  in  its  application? 

8.  Are  the  assumed  stress  or  displacement  functions  simple? 

9.  Is  the  derivation  of  the  element  matrices  simple,  straightforward  and 
understandable,  and  are  the  matrices  given  in  explicit  form? 

10.  What  element  matrices  are  available? 

11.  Is  the  element  a  synthesized  element,  that  is,  is  it  built-up  of  a 
number  of  elements? 
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12.  Are  the  element  matrices  easy  to  code? 

13.  Is  the  element  theory  and  coding  well-documented  together  with  check 
examples? 

14.  Does  the  element  contain  any  false  zero  energy  modes  (false  zero  stress 
or  strain  states)  other  than  those  corresponding  to  rigid  body  modes? 

15.  Is  the  element  mixable  with  other  elements  within  a  particular  program? 

16.  Can  the  element  be  easily  incorporated  into  a  particular  program? 

17.  Has  the  element  been  fully  checked  out?  How  often  and  in  what  problems 
has  it  been  used?  The  debugging  of  an  element  is  very  difficult  and 

is  made  even  more  difficult  by  element  complexity. 

18.  Is  the  originator  of  the  element  available  for  consultation  if  problems 
arise?  In  many  cases,  the  software  developers  are  unaware  of  the 
detailed  theory  behind  their  elements. 

19.  What  are  the  element  material  properties? 

20.  Does  the  element  help  to  reduce  the  effort  to  prepare  the  model,  to 
obtain  computer  results,  to  interpret  the  results? 

21.  For  a  given  model,  do  the  results  change  if  the  element  specifying 
nodes  are  written  in  a  different  order? 

Table  7-1  is  our  assessment  of  how  COSMIC/NASTRAN  rates  in  these  categories  for 
the  selected  group  of  elements.  This  assessment  is,  in  many  cases,  quite 
subjective  and  is  based  largely  on  the  documentation  (primarily  the  Theoretical 
Manual).  Results  show  that  in  general,  according  to  Robinson's  criteria, 
COSMIC/NASTRAN  rates  a  "fair"  ranking  in  its  documentation  for  any  particular 
element. 

Robinson  also  recommended  these  basic  evaluation  tests  (References  31,34)  ; 

•Single  element  test 

•  Convergence  test 

»  Element  completeness  test 

•  Patch  test 

•  False  zero  energy  mode  test 


I 
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The  "single  element  test"  consists  of  taking  a  single  element  in 
rectangular  form,  considering  it  as  a  structure  and  then  investigating  its 
behavior  under  various  loading  and  constraint  conditions  for  various  geometries. 
This  test  is  very  helpful  in  understanding  the  load-carrying  capabilities  of 
the  element  and  shows  its  sensitivity  to  aspect  ratio.  (This  test  is  relatively 
simple  to  apply  only  for  triangular  or  quadrilateral  membrane  and  plate 
bending  elements) . 

The  "convergence  test"  consists  of  studying  the  behavior  of  some  parameter 
or  parameters  for  various  numbers  of  elements  in  a  model  for  a  particular  problem. 
As  the  number  of  elements  in  a  model  is  increased,  and  the  mesh  becomes  more 
refined,  the  results  should  approach  the  true  solution. 

The  "element  completeness  tests"  consist  of  showing  that  the  element  dis¬ 
placement/stress  assumptions  contain  the  constant  strain/stress  modes  and  the 
zero  strain/stress  (rigid  body)  modes.  These  are  essential  requirements  for 
any  finite  element.  If  each  element  in  a  model  passes  these  tests,  and  the  model 
is  a  compatible/equilibrium  one,  the  results  will  converge  monotonically. 

The  "patch  test"  consists  of  applying  a  constant  strain/stress  state  to 
an  assemblage  of  elements  and  showing  that  all  the  elements  contain  the  same 
constant  strains/stresses.  This  test  is  for  "nonconforming"  or  "noncompatible" 
and  "non-equilibrium"  models  and  should  be  passed  if  the  results  are  to  converge 
(though  not  necessarily  monotonically).  The  patch  test  is  a  completeness  con¬ 
dition  for  an  assemblage  of  elements  as  against  an  individual  element.  The 
primary  purpose  of  this  test  is  to  provide  an  estimate  of  the  error  which  will 
be  incurred.  It  is  noted  that  increasing  the  number  of  elements  in  the  "patch" 
will  probably  reduce  the  error  if  the  element  performs  well  in  the  convergence 
test. 
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The  "false  zero  energy  mode  test”  consists  of  showing  that  an  element  does 
not  contain  zero  strain/stress  fields  other  than  those  associated  with  rigid 
body  modes.  Elements  can  contain  false  zero  states  which  are  due  to  the  basic 
strain/stress  assumptions  and  the  subsequent  choice  of  nodal  degrees  of  freedom. 

In  addition  to  these  tests,  we  will  often  include  the  results  on  "element 
efficiency  test".  This  test  will  in  general  be  different  for  each  type  of 
element,  and  shows  a  measure  of  element  performance  and  coding  efficiency  for 
that  problem.  Some  well-known  benchmark  cases  used  by  many  researchers  include: 

.Plate  bending  element.  Simply  supported,  or  clamped,  square  plate  loaded 
centrally  by  a  concentrated  load. 

.Plate/shell  element.  Scordelis-Lo  cylindrical  shell  problem  loaded  by  its 
own  weight. 

.Conical  shell  element.  Cantilevered  cylindrical  shell  loaded  at  its  free 
end  by  a  line  load  around  the  circumference. 

.Solid  of  revolution  element.  Thick-walled  cylinder  loaded  by  internal 
pressure. 

.3-D  solid  elements.  Slender  cantilever  beam  problem,  a  pressurized  thick- 

walled  cylinder  problem,  and  a  study  of  element  eigenvalues  and  eigenvectors. 

At  the  end  of  each  element  evaluation,  an  assessment  will  be  made  of  that  element's 
quality  and  efficiency. 
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7.2.  The  QDPLT  Element 

QDPLT  was  selected  for  advanced  evaluation  because  it  is  one  of  the  more 
complex  plate  bending  elements  in  the  NASTRAN  library.  Also,  there  is  less 
documentation  in  the  Theoretical  Manual  for  this  element  than  for  most  other 
shell  elements.  In  the  course  of  this  investigation,  it  was  determined  that 
deficiencies  had  been  previously  observed  in  QDPLT  (Ref.  33) .  QDPLT  consists 
of  four  nonconforming  TRPLT's,  models  bending  behavior  only,  and  does  not  have 
membrane-bending  coupling. 

7.2.1.  Single  Element  Test 

The  single  element  test  for  a  plate  bending  element  consisted  of  evaluating 
the  displacement  and  rotation  for  eight  cases  (giving  various  combinations  of 
constraints  and  loads)  for  a  range  of  element  aspect  ratios: 

Cases  1,  2,  3,  and  4:  Displacement  in  the  z-direction  (U  )  and  rotation 

z 


Cases  5,6,7,  and  8: 


about  the  y-axis  at  node  3  (0^) . 

Displacement  in  the  z-direction  (Uz)  and  rotation 
about  the  x-axis  at  node  3  (0^) . 


The  single  element  test  results  are  plotted  in  Figures  7.2-1  to  7.2-8,  where 
they  are  compared  with  the  standard  beam  solutions. 

For  our  evaluation  of  QDPLT,  the  following  material  properties  were  used: 
Young's  modulus  E  =  107  psi 


Shear  modulus 


2 (1+v) 


=  3.8462  x  10  psi 


Poisson's  ratio  v  =  0.3 


For  example,  the  theoretical  values  are  presented  on  the  following  three  pages: 


CASE  BOUNDARY  CONDITIONS  THEORETICAL  RESULTS 


CASE  .  BOUNDARY  CONDITIONS  THEORETICAL  RESULTS 
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The  critical  single  element  tests  for  QDPLT  are  Cases  6  and  8  (Figures  7.2-6 
and  -8).  In  both  cases  where  the  performance  was  poor,  a  twisting  moment  was 
applied  to  the  element  with  one  edge  fully  clamped.  This  activates  differential 
(bilinear)  bending.  It  is  interesting  to  note  that  when  the  rotational  restraint 
was  not  applied  (Cases  5  and  7)  with  the  loading  otherwise  identical,  the  element 
performed  well. 

Other  4-node  plate  bending  elements,  such  as  Lockheed/NASTRAN' s  LORA  and  MSC/ 
NASTRAN's  QUAD 4  correlated  reasonably  well  with  theory  up  to  aspect  ratios  of 
12  (Ref.  33).  Except  for  load  case..  6  and  8, QDPLT  performed  well  in  the  single 
element  tests. 


7.2.2.  Convergence  Test 

A  simple  convergence  test  often  used  for  plate  bending  elements  is  a  square 
plate  centrally  loaded  by  a  concentrated  load  P,  with  simply- supported  and 


Theoretical  deflections  at  the 

center  are  (Ref.  36) : 

aPb2 

6  = 

3 

Assume: 

b  =  2  in. 

Et 

P  =  4  lb. 

^2. 0272  x  10-4  in. 

s imply- s uppor ted 

t  =  0.1  in. 

6  =  < 

1  -5 

E  =  107  psi 

[9.976  x  10  in. 

clamped 

mm 


FIGURE  7.2-1.  SINGLE  ELEMENT  TEST  FOR  QDPLT.  CASE  1.  FIGURE  7.2-2.  SINGLE  ELEMENT  TEST  FOR  QDPLT.  CASE  2 


COMPUTED  CENTRAL  DEFLECTION  (in.) 


NO. 

OF  ELEMENTS 

SIMPLY 

SUPPORTED 

CLAMPED 

FOR 

1/4  OF  PLATE 

HBHRH 

%  Error 

5  x  105  (in.) 

%  Error 

1 

2.198274 

8.43 

5.475623 

-42.3 

4 

2.134809 

5.31 

9.989985 

0.140 

9 

2.109106 

4.04 

10.33296 

3.57 

16 

2.100579 

3.61 

10.41435 

4.39 

These  QDPLT  results  are  plotted  in  Figures  7. 2. 2-1  and  7. 2. 2-2,  and  show  the  same 
convergence  trends  as  in  Ref.  31.  Unfortunately,  the  modeling  error  discussion 
(Sec.  15.2)  in  Ref.  7  does  not  contain  QDPLT  results.  It  is  interesting  to  note 
that  the  convergency  trend  is  not  monotonic  and,  in  fact,  appears  to  diverge 
slightly  for  the  clamped  case  as  the  number  of  elements  is  increased,  with  four 
elements  yielding  very  accurate  results. 


CENTRAL  DEFLECTION  (INCHES 


2.2 


2.06, 

2.04 

2.02 


THEORY 


NUMBER  OF  ELEMENTS 

FIGURE  7. 2.2-1.  CONVERGENCE  TEST  FOR  QDPLT.  SIMPLY  SUPPORTED  SQUARE  PLATE 
(AR  -  1)  WITH  CENTRAL  POINT  LOAD. 


FIGURE  7. 2. 2-2.  CONVERGENCE  TEST  FOR  QDPLT.  CLAMPED  SQUARE  PLATE  (AR  =  1) 


WITH  CENTRAL  POINT  LOAD 
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7.2.3.  Patch  Test 


The  patch  test,  first  proposed  by  B.  Irons  in  1972,  requires  a  nonconforming 
element  to  produce  constant  stresses/strains  when  these  elements  are  formed  to¬ 
gether  in  a  model  which  has  at  least  one  internal  node.  It  thus  measures  completeness 
for  an  assemblage  of  elements.  For  a  quadrilateral  plate  bending  element  such  as 
QDPLT ,  one  possible  patch  test  is  this  “picture  frame"  model  subjected  to  the 
loads  shown  (Ref.  31): 


When  we  used  such  a  model  for  QDPLT,  and  assuming: 

7 

a  =  b  =  0.5,  E  =  10  psi,  v  =  0.3,  t  =  0.1  in. 
We  obtained  stress  results  as  follows: 

Elements  1,  2,  4,  5:  a  =  1421  psi 


Element  3: 


a  =  1138  psi 
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QDPLT  therefore  fails  the  patch  test.  The  stress  values  given  are  the  Von 
Mises  effective  stresses  at  the  centroid  of  the  element.  The  Von  Mises 
effective  stress  was  chosen  as  a  measure  of  stress  level  because  the  orien¬ 
tation  of  the  element  i-j  nodes  defines  the  stress  output  orientation  and 
therefore  stress  components  are  not  directly  comparable  from  element  to 
element.  This  result  indicates  that  an  assembly  of  elements  will  not  pro¬ 
duce  identical  stresses  in  each  element  even  though  the  applied  loads  should 
produce  a  uniform  stress  rate.  This  is  not  an  unexpected  result,  however, 
since  slope  compatibility  is  not  satisfied  at  the  mutual  boundaries  of  the 
elements . 

In  order  to  determine  which  type  of  loading  produced  the  most  severe 

disparity,  the  combined  loading  was  divided  into  individual  load  components 

as  shown  in  Figure  7. 2. 3-1.  As  can  be  seen,  the  loading  consists  of  edge 

moments,  M  and  M  ,  and  an  out-of-plane  force  at  one  corner.  The  two  moment 
x  y 

load  cases  should  produce  a  constant  curvature,  R  or  R  ,  resulting  in  a 

x  y 

uniform  stress  distribution  on  each  surface  of  the  plate.  The  third  loading 
condition  produces  a  twisting  load  on  the  plate  and  should  result  in  a  constant 
state  of  stress  in  all  elements.  The  results  indicate  that  the  moment  loads 
produce  the  largest  range  of  stress  values  in  the  elements. 

A  detailed  examination  of  the  displacement  solution  revealed  that  the 
slopes  were  discontinuous  across  the  element  boundaries,  therefore  the  constant 
curvatures  were  not  correctly  predicted.  It  was  interesting  to  note  that  the 
concentrated  load  at  the  corner  resulted  in  a  fairly  uniform  stress  field, 
even  though  this  appears  to  be  the  most  severe  loading  condition  due  to  the 
warping  of  the  element.  However,  the  single  element  test,  Figure  7.2-7,  does 
indicate  that  the  element  will  perform  well  under  this  loading  as  long  as 
moment  constraints  are  not  applied  (as  in  Case  8,  Figure  7.2-8).  Possibly, 
the  reason  for  the  poor  performance  of  the  QDPLT  element  in  modeling  the  con¬ 
stant  curvature  is  related  to  its  poor  performance  in  Cases  6  and  8  of  the 
single  element  tests,  where  slope  constraints  were  enforced  at  the  support. 

The  poor  performance  of  QDPLT  in  the  patch  test  illustrates  the  need  for 
the  analyst  to  be  aware  of  the  limitations  of  QDPLT  or  any  other  plate  bending 
elements  that  are  being  used  in  a  structural  analysis.  The  results  should 


I 
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be  checked  thoroughly,  especially  in  areas  where  large  stress  or  strain 
gradients  occur.  Also,  the  fineness  of  the  mesh  should  be  such  that  the 
difference  in  stress  between  adjacent  elements  is  relatively  small.  Areas 
of  the  mesh  that  exhibit  large  stress  differences  between  elements  should 
either  be  reanalyzed  with  a  finer  mesh  or  results  should  be  assigned  a 
fairly  generous  error  range  when  the  stresses  are  used  to  assess  the  adequacy 
of  the  structure. 


f 
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N  -  Node 
©  -  Element 


A\x  =  .5 


INDIVIDUAL  LOADING  CONDITION 


CONSTANT  STRESS  SHAPES 


ALL  LOADS  APPLIED 
SIMULTANEOUSLY 


z 


VON  MISES  EFFECTIVE  STRESS  AT  CENTROID 


FIGURE  7. 2. 3-1.  QDPLT  PATCH  TEST  RESULTS  FOR  INDIVIDUAL  LOADING  CONDITIONS. 
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7.2.4.  Element  Completeness  Test 

The  completeness  requirement  for  an  element  is  satisfied  if  the  rigid  body 
displacements  and  constant  strain  states  are  represented.  A  rigid  body  displace¬ 
ment  is  the  most  elementary  deformation  that  an  element  may  undergo.  This 
condition  states  that  there  should  exist  combinations  of  values  of  the  generalized 
coordinates  that  cause  all  points  on  the  element  to  experience  the  same 
displacement.  One  such  combination  should  occur  for  each  of  the  rigid  body 
translations  and  rotations.  In  a  typical  displacement  model,  the  constant  term 
q^  provides  for  a  rigid  body  displacement. 

The  constant  strain  requirement  can  be  stated  in  similar  terms.  There  should 
exist  combinations  of  values  of  the  generalized  coordinates  q^  that  cause  all 
points  on  the  element  to  experience  the  same  strain.  As  a  body  is  subdivided 
into  smaller  and  smaller  elements,  the  strains  in  each  infinitesimal  element 
approach  constant  values.  Unless  these  constant  strains  are  included  in  the  shape 
function,  convergence  to  the  correct  solution  will  not  occur.  In  a  typical 
displacement  model,  the  terms  (q2*)  and  q^y)  provide  for  constant  strains 
and  £y  in  an  element. 

QDPLT  is  a  quadrilateral  bending  element  which  is  composed  of  four  overlapping 
basic  bending  triangles.  Upon  close  examination  of  the  theoretical  development 
of  this  element  in  Ref.  7,  it  was  concluded  that  the  terms  for  both  rigid  body 
displacements  and  constant  strain  states  are  included  in  the  displacement  model. 

7.2.5.  False  Zero  Energy  Mode  Test 

If  the  number  of  element  independent  deformation  variables  plus  the  number 
of  rigid  body  modes  is  less  than  the  number  of  nodal  degrees  of  freedom  in  the 
local  system,  an  element  may  contain  falze  zero  energy  modes  (Ref.  31,34).  Zero 
energy  modes  are  always  associated  with  rigid  body  modes  since  these  give  zero 
elastic  energy  and  zero  stress  or  strain  fields.  Even  when  the  number  of  inde¬ 
pendent  deformation  variables  plus  rigid  body  modes  is  equal  to  the  number  of 
nodal  degrees  of  freedom,  false  modes  can  exist. 

False  zero  energy  modes  can  be  found  by:  numerical  techniques  directly  using 
the  element  matrices;  an  iterative  static  method  of  gradually  increasing  the 
number  of  support  freedoms  (Ref.  31);  or,  a  modal  procedure  whereby  the  eigen¬ 
values  are  computed  for  an  unsupported  structure  (Ref.  34) ,  as  suggested  by 
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Wilson,  Taylor,  and  Doherty  (1969).  Although  a  single  element  may  contain 
false  zero  energy  modes,  in  many  cases  experience  has  shown  that  these  modes 
do  not  propagate  when  combined  with  other  elements  to  form  a  structural  model, 
therefore  causing  no  problems.  To  perform  this  test  rigorously  for  QDPLT  would 
have  meant  setting  up  either:  a  static  method  of  iteratively  searching  for  false 

modes  by  gradually  increasing  the  number  of  support  freedoms - for  a  variety  of 

supports  and  loading  cases;  or  an  exhaustive  eigenvalue  interpretation  procedure 
similar  to  what  was  done  for  the  NASTRAN  solid  elements  (Subsection  7. 7. 2. 3). 
Because  of  the  relative  lack  of  importance  of  this  test  and  our  extensive  work 
already  previously  reported  to  characterize  QDPLT  performance  in  the  single 
element  test,  convergence  test,  and  patch  test,  it  was  decided  to  forego  the 
false  zero  energy  mode  test  on  QDPLT. 

7.2.6.  Conclusions 

Advanced  evaluation  tests  on  the  QDPLT  bending  element  show  it  rated  fair  in 
performance.  In  summary,  it  may  be  said  that  QDPLT:  performed  well  in  the  single 
element  test  in  six  out  of  eight  cases  of  various  loading/support  combinations 
(the  exceptions  being  the  twist  load  cases  6  and  8);  converged  satisfactorily  in 
the  two  convergence  tests  of  a  simply- supported  and  clamped  square  plate  under  a 
concentrated  center  load;  and  failed  the  patch  test.  Its  overall  performance 
is  inferior  to  other  currently  available  4-noded  plate  elements,  such  as: 
MSC/NASTRAN ' s  QUAD4 ,  Lockheed-California  NASTRAN' s  LORA  stress-based  plate  bending 
element,  and  T.  H.  H.  Pian's  stress-based  element  with  nine  independent  force 
variables  (Ref.  33). 
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7.3.  The  TRSHL  Element 

TRSHL  is  a  triangular  shallow  shell  element  implemented  into  Level  16.0 
of  COSMIC /NAS TRAN  in  1976.  When  TRSHL  was  originally  selected  for  advanced 
evaluation,  it  seemed  a  logical  choice  because  it  is  a  higher-order  shell 
element  recently  incorporated  into  the  COSMIC/NASTRAN  element  library.  In 
this  subsection,  a  description  of  TRSHL  documentation,  theory,  status, 
and  efficiency  will  be  given.  Most,  or  all  of  this  background 
information  is  unavailable  to  the  average  NASTRAN  user,  and  it  is  included  to 
illustrate  the  implementation  history  and  current  status  of  a  typical  element 
in  the  NASTRAN  library. 

7.3.1.  Background 

TRSHL  was  developed  and  implemented  into  COSMIC/NASTRAN  by  Dr.  R.  Nara- 
yanaswami  in  1974-1976  (References  37,38)  while  he  was  at  Old  Dominion 
University,  Norfolk,  Virginia.  The  element  models  curved  thin  shell  behavior 
by  approximating:  the  membrane  behavior  using  the  linear  strain  TRIM6  element, 
the  bending  behavior  using  the  quintic  transverse  displacement  TRPLT1  element, 
and  membrane-bending  coupling  using  Novozhilov  shallow-shell  theory.  This 
6-noded  element  has  30  degrees  of  freedom;  at  each  node  there  are  five  degrees 
of  freedom  consisting  of  three  translations  and  two  rotations.  The  element 
is  designed  for  use  with  the  statics,  normal  modes,  and  buckling  rigid  formats 
of  NASTRAN. 

In  our  early  investigations,  two  interesting  things  were  discovered  about 
TRSHL.  First,  the  documentation  of  TRSHL  in  the  NASTRAN  manuals  was  poorly 
written,  short,  and  incomplete.  Narayanaswami ' s  two  Old  Dominion  University 
reports  cited  in  the  NASTRAN  Theoretical  Manual  were  unavailable  through  the 
usual  NASA  or  NTIS  publication  sources,  even  though  the  research  was  funded 
by  NASA/Langley.  (We  were  able  to  obtain  Reference  37  by  writing  directly 
to  Old  Dominion  University).  Second,  although  TRSHL  still  exists  in  the 
COSMIC/NASTRAN  Level  17.5  element  library,  its  use  appears  to  be  limited  and 
has  been  abandoned  in  both  the  latest  versions  of  NASTRAN  offered  by  the 
MacNeal-Schwendler  Corporation  (MSC)  and  Universal  Analytics,  Inc.  (UAI) . 

The  history  of  the  development  and  implementation  of  the  TRSHL  element 
illustrates  one  of  the  weaknesses  of  the  piecemeal  development  of  a  large 
general  purpose  program  such  as  NASTRAN.  The  use  of  a  number  of  different 
subcontractors  to  provide  finite  elements  or  solution  algorithms  often  leads 
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to  a  situation  where  the  developer  is  not  available  to  completely  debug 
the  program  modification  in  its  final  form.  Furthermore,  the  complete¬ 
ness  and  accuracy  of  the  documentation  often  vary  significantly. 

Additionally,  the  developer  is  usually  not  available  to  provide  user 
support  for  this  part  of  the  program.  This  latter  criticism  is  often 
cited  (for  example  Ref.  22)  as  an  important  evaluation  point  in  comparing 
computer  programs. 

Universal  Analytics,  Inc.  was  contacted  to  inquire  why  TRSHL  is 
not  included  in  UAI/NASTRAN  anymore.  In  response,  Mr.  M.  J.  Morgan  of  UAI 
informed  us  that  UAI  found  TRSHL  did  have  poor  performance  and  reputation, 
and  UAI  had  decided  instead  to  implement  its  own  QUAD3  and  QUAD8  elements 
for  shell  analysis.  He  referred  us  to  Dr.  A.  B.  Potvin  of  Exxon  Production 
Research  Company  (Houston,  Texas)  for  more  TRSHL  evaluation  results.  Dr.  Potvin 
told  us  that  he  discovered  several  minor  coding  errors  in  the  TRSHL  implemen¬ 
tation  into  Level  16.0  of  COSMIC/NASTRAN ,  and  after  correcting  these  himself, 
he  found  its  performance  and  convergence  to  be  relatively  poor  in  the 
Scordelis-Lo  problem  (Ref.  39).  Potvin’s  TRSHL  data  are  included  in  Figure 
7.3-2. 

Based  upon  these  inquiries,  the  following  conclusions  have  been  drawn 
concerning  TRSHL:  The  documentation  in  NASTRAN  is  inadequate;  its  use  is 
apparently  quite  limited;  there  is  little  technical  support  available  from  its 
developer;  its  performance  is  inferior  (as  can  be  seen  from  Figure  7.3-2) 
to  other  currently  available  shell  elements,  and  its  coding  implementation 
into  COSMIC/NASTRAN  was  apparently  not  thoroughly  tested. 

7.3.2.  Theory 

The  TRSHL  element  uses  only  displacements  and  rotations  as  nodal  degrees 
of  freedom  (Figure  7.3-1).  The  displacement  functions  for  u,  v,  and  w  are 
chosen  as  follows: 

2  2 

u  =  qj  +  q2x  +  q3y  +  q4x  +  q5xy  +  qgy 

V  =  q7  +  q8x  t  q9y  +  q^  +  quxy  ♦  q^ 


FIGURE  7.3-1.  TRSHL  SHELL  ELEMENT  GEOMETRY  AND  COORDINATE  SYSTEMS 


7-28 


The  element  includes  transverse  shear  flexibility  in  the  stiffness  formulation. 
Use  of  quadratic  polynomials  for  the  shell  surface  geometry  implies  constant 
curvatures,  thereby  agreeing  with  the  approximations  of  shallow  shell  theory. 

7.3.3.  Evaluation  Results 

Although  Narayanaswami  presented  two  numerical  examples  (spherical  cap, 

Scordelis-Lo  cylindrical  shell),  here  we  will  comment  on  only  the  latter.  Table 

7.3-1  shows  TRSHL  results  for  the  Scordelis-Lo  problem,  and  the  trends  of  four 

predicted  displacements  and  three  stress  resultants  versus  mesh  size.  Note  that 

even  for  the  fine  6x6  mesh,  discrepancies  are  quite  pronounced  for  four  of  the 

seven  results:  V„,  W„,  M  ,  M 

B  C  yy  xx 
JJc  c 

Narayanaswami  (Ref.  37)  did  not  provide  a  complete  explanation  for  the 
difference  between  the  calculated  values  and  the  exact  solution.  If  one  examines 
only  the  free  edge  midpoint  vertical  displacement  (Figure  7.3-2),  TRSHL  fared 
somewhat  better,  though  it  is  still  inferior  to  Cowper,  et.  al.  (1970)  and  many 
other  available  commercial  shell  elements.  Potvin's  independent  check  of  TRSHL 
results  (Ref.  39)  is  also  shown:  it  is  noted  that  there  is  a  considerable  dif¬ 
ference  between  Potvin’s  results  and  those  reported  by  Narayanaswami.  It  is  also 
interesting  to  note  that,  although  TRSHL  is  described  in  Reference  7,  it  was 
omitted  in  a  table  which  lists  available  COSMIC/NASTRAN  finite  elements  in 
Reference  11. 

7.3.4.  Conclusions 

Based  on  the  results  of  this  investigation,  it  is  concluded  that  TRSHL 
element  quality  and  efficiency  rate  considerably  below  those  of  currently 
available  shell  elements  in  competing  software.  A  low  (and  declining)  level 
of  its  use  in  the  future  is  predicted,  leading  to  its  probable  ultimate  demise 
in  the  COSMIC/NASTRAN  element  library. 


TABLE  7.3.1.  TRSHL  RESULTS  FOR  SCORDELIS-LO  CYLINDRICAL  SHELL  ROOF  PROBLEM 
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7.4.  The  TRIM6  and  TRPLT1  Elements 

Closely  related  to  TRSHL  are  the  linear  strain  triangular  membrane  element 
TRIM6  and  the  higher  order  triangular  bending  element  TRPLT1.  This  section  gives 
a  short  summary  of  the  element  theory  and  efficiency  for  both  elements.  Evaluation 
results  are  already  well  documented  in  References  7  and  38.  TRIM6,  TRPLT1,  and 
TRSHL  were  added  to  Level  16.0  of  COSMIC/NASTRAN  by  Dr.  Narayanaswami  in  December 
1976. 

7.4.1.  Theory 

The  linear  strain  triangular  membrane  element  TRIM6  is  based  on  J.  H.  Argyris* 
derivation  in  1965.  The  element  has  six  nodes,  three  at  the  vertices  and  three 
at  the  midpoints  of  the  sides.  The  element  uses  a  quadratic  displacement  field. 

2  2 

u  =  +  q2x  +  q3y  +  q4x  +  q5xy  +  qgy 

v  =  q?  +  q0x  +  q9y  ♦  q^X2  +  q^y  +  q^y2 

The  thickness  of  the  element,  as  well  as  the  temperature  distribution  within  the 
element,  can  have  bilinear  variation.  TRIM6  may  be  used  for  the  statics  and  normal 
modes  rigid  formats,  but  cannot  be  used  for  differential  stiffness  and  piecewise 
linear  analyses.  Element  stresses  are  computed  at  the  three  vertices  and  at  the 
centroid. 

TRPLT1  was  developed  by  Narayanaswami  in  1974  as  a  modification  of  Cowper  et.  al. 
(1968) .  Like  TRIM6  and  TRSHL,  TRPLT1  has  six  nodes,  three  at  the  vertices  and  three 

at  the  midpoints  of  the  sides.  A  quintic  displacement  field  is  chosen  for  the 

transverse  displacement  w  (see  Subsection  7.3.2).  The  element  is  nonconforming, 
meaning  it  has  no  slope  continuity  for  two  adjacent  elements  with  a  common  edge. 

Each  node  has  three  degrees  of  freedom:  two  rotations  and  the  transverse  displace¬ 
ment  w.  Transverse  shear  flexibility  is  included  in  the  stiffness  formulation. 

The  element  thickness  can  have  bilinear  variations.  TRPLT1  may  be  used  in  the 
statics  and  normal  modes  rigid  formats.  The  element  internal  forces  are  recovered 
at  the  three  vertices  and  at  the  centroid. 

7.4.2.  Evaluation  Results 

The  performances  of  TRIM6  and  TRPLT1  are  already  adequately  described  in 
References  7  and  38.  We  decided  not  to  repeat  these  exercises,  and  only  report 
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here  some  representative  test  problems  which  demonstrate  their  efficiency  and 
convergence  characteristics. 

.Example  1.  Cantilevered  beam  using  membrane  elements. 

Figure  7.4-1  shows  the  excellent  correlation  of  TRIM6  with  theory  and  its 
superiority  over  the  QDMEM1  isoparametric  membrane  element  and  the  conventional 
quadrilateral  membrane  element  QDMEM.  This  example  demonstrates  the  improved 
accuracy  due  to  use  of  a  quadratic  displacement  polynomial,  and  appears  in 
both  References  7  and  38. 

.Example  2.  Rectangular  plate  bending  using  TRPLT1. 

Figure  7.4-2  gives  a  typical  TRPLT1  deflection  convergence  study  for  a 
rectangular  plate,  simply  supported  at  the  four  edges  and  with  a  central 
concentrated  load.  Note  that  TRPLT1  converges  monotonically  from  above  the 
theoretical  value,  a  behavior  which  is  typical  of  nonconforming  elements 
and  suggesting  need  for  a  patch  test.  This  example  was  for  a  (b/a)  ratio  of  2. 

7.4.3.  Conclusions 

A  careful  review  of  all  the  convergence  efficiency,  and  demonstration  examples 
in  References  7  and  38  reveals  that  the  TRIM6  and  TRPLT1  elements  performed  quite 
well  for  a  variety  of  conditions  and  load  cases  -  including  statics,  thermal  loads, 
vibration,  and  buckling  problems. 
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FIGURE  7.4-2.  CENTRAL  DEFLECTION  OF  RECTANGULAR  PLATE 


7.5.  The  CONEAX  Element 


CONEAX  is  a  conical  shell  element  which  can  have  either  axisymmetric  or 
asymmetric  loading.  This  section  briefly  explains  CONEAX  theory  and  our 
advanced  evaluation  results  using  an  edge-loaded  cylinder  to  assess  its  accuracy. 

7.5.1.  Theory 

The  CONEAX  element  coordinate  geometry  and  5  degrees  of  freedom  are  shown 
below. 


The  displacements  u  and  v  are  assumed  to  be  linear  functions  of  the  meridional 
position  s,  while  w  is  assumed  to  be  a  cubic.  For  each  harmonic  order  q^,  the 
harmonic  components  of  deflection  are  {Ref.  7): 


VS)  =  qln 


+  q2nS 


v  (s)  =  q_  +  q.  s 
n  3n  Mn 


2  3 

w  (s)  =  q_  +  q,  s  +  q„  s  +  q„  s 
n  7  ^5n  ^6n  ^7n  M8n 


Therefore,  CONEAX  is  similar  in  derivation  to  the  conical  shell  element  of  Grafton 
and  Strome  (Ref.  44) .  In  a  conical  assemblage  the  structure  modeled  by  this  conical 
element  is  permitted  to  deform  only  approximately  but  is  loaded  in  a  manner  con¬ 
sistent  with  those  deformations.  The  philosophy  is  that  many  small  elements,  each 


j 


[ 

deforming  in  a  simple  way,  should  provide  a  good  approximation  to  the  actual 
structure.  In  any  region  of  the  structure  in  which  the  character  of  the  deformation 
varies  rapidly,  it  is  necesrary  to  use  a  fairly  large  number  of  such  elements. 
Clearly,  this  is  the  case  in  edge  bending  calculations  (Ref.  45). 

Properties  of  CONEAX  are  assumed  axisymmetr ic.  However,  the  loads  and  de¬ 
flections  need  not  be  axisymmetric;  they  are  expanded  in  Fourier  series  with 
respect  to  the  azimuth  coordinate.  An  unusual  feature  of  CONEAX  is  that  it  includes 
transverse  shear  flexibility.  The  rotations  a  and  g  are  independent  motions 
because  of  the  transverse  shear  flexibility.  Rotation  about  the  normal  to  the 
surface  is  not  included,  since  such  rotation  can  be  adequately  represented  by 
the  gradients  of  u  and  v.  The  element  cannot  be  combined  with  other  NASTRAN 
structural  elements  in  the  solution  of  problems. 

7.5.2.  Advanced  Evaluation 

To  evaluate  the  convergence  and  accuracy  of  CONEAX,  the  radial-loaded  cylinder 
example  used  in  Grafton  and  Strome's  paper  (Ref.  44)  and  by  many  subsequent  re¬ 
searchers  was  selected  as  the  benchmark.  This  example  features  rapidly  varying 
deflections  and  stress  resultants. 
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Theoretical  predictions  for  the  radial  deflection  and  meridional  moment  can  be 
obtained  using  Hetenyi  (Ref.  46) : 

0o  -fix 

Radial  deflection:  w  =  — —  e  cos 

28  D 


Meridional  moment: 


Qo  -gx  . 

M  =  -—  e  sin  Bx 

8 


where : 
( V= • 3 ) 


1.285 

/Rt 


1.285 

/5(.01) 


5.746694 


D 


Et“ 


12 (1-v  ) 


7  3 

10  (.01) 

12(1-. 32) 


0.915751 


Computed  results  using  CONEAX  were  remarkably  close  to  theory,  even  for  Case  1 
(5  elements) .  All  three  models  produced  consistently  better  correlation  with 
theory  for  radial  deflection  and  meridional  moment,  as  compared  to  Grafton  and 
Strome  (Ref.  44).  Results  are  listed  below  and  plotted  in  Figures  7.5-1  and 
7.5-2. 
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.Meridional  Moment:  (x  is  at  element  centroid) 
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FIGURE  7.5-1.  RADIAL  DEFLECTION  -  CONEAX  VS.  THEORY 
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7.5.3,  Other  CONEAX  Uses 

In  addition  to  handling  axisymmetric  loadings,  CONEAX  is  versatile  in 
its  use  and  can  be  applied  to  a  variety  of  other  loading  conditions: 

•Unsymmetric  loading  -  using  standard  Fourier  expansion  techniques 

•Line  loads  -  uniformly  distributed  or  harmonic  coefficients 

•Concentrated  loads  -  at  specified  azimuth  positions 

•Pressure  loads  -  normal  to  surface 

•Gravity  loads 

•Thermal  loads 

•Enforced  strains 

•Enforced  displacements 

♦Branched  shells 

7.5.4.  Conclusions 

Our  investigations  indicate  that  CONEAX  is  a  good  conical  shell  element, 
with  good  convergence  properties  and  exceptional  versatility.  Only  the 
axisymmetric  mechanical  load  capability  was  evaluated  in  this  study.  It  has 
been  reported  (Ref.  53)  that  the  thermal  loading  capability  that  exists  may 
contain  errors  and  is  not  recommended  for  use. 


I 
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7.6.  The  TRAPRG  and  TRAPAX  Elements 

TRAPRG  is  one  of  the  four  solid-of-revolution  elements  offered  in  COSMIC/ 
NASTRAN.  It  is  a  four-noded  axisymmetric  trapezoidal  ring  element  with  axi- 
symmetric  loading  capability.  The  other  three  are:  TRIARG,  same  as  TRAPRG  except 
it  is  triangular;  TRAPAX  and  TRIAAX,  axisymmetric  trapezoidal  and  triangular 
ring  elements  with  nonaxisymmetr ic  loading  capability.  The  benchmark  problem 
selected  to  test  TRAPRG  and  TRAPAX  performance  is  a  thick-walled  cylinder  under 
internal  pressure,  one  used  by  many  past  researchers,  and  a  good  test  to  measure 
their  capability  to  represent  simple  radial  dilatation  behavior. 

7.6.1.  Theory 

According  to  the  Theoretical  Manual  (Ref.  7)  ,  the  formulation  of  the  theory 
behind  TRAPRG  is  "mathematically  consistent"  with  that  in  Mallet  and  Jordan's 
MAGIC  (1969)  and  MAGIC  II  (1971)  codes  developed  at  the  Air  Force  Flight  Dynamics 
Laboratory.  The  two  original  finite  element  papers  on  the  subject  of  structural 
analysis  of  axisymmetric  solids  were  by  Clough  and  Rashid  (1965)  and  Wilson  (1965) 
There  was  apparently  (Ref.  7)  an  extension  of  the  Wilson  formulation  to  include 
orthotropic  material  properties  for  TRAPRG  in  COSMIC/NASTRAN .  The  TRAPRG  theo¬ 
retical  derivations  are  quite  standard  for  solid-of-revolution  elements.  The 
31-page  documentation  in  Reference  7  is  adequate  in  describing:  coordinate  no¬ 
tations;  displacement  functions  (Fourier  harmonics  in  circumferential  direction); 
derivation  of  the  element  stiffness  matrix  and  mass  matrix;  treatment  of  ortho¬ 
tropic  behavior;  development  of  non-axisymmetr ic  load  vectors  (including  thermal 
loads);  prestrain,  pressure  and  thermal  load  vectors;  displacement  and  stress 
recovery.  The  TRAPRG  element  and  the  other  axisymmetric  elements  are  therefore 
ten  to  fifteen  years  old.  They  represent  fairly  typical  solid-of-revolution 
elements  available  in  all  general  purpose  codes  and  documented  in  every  finite 
element  text,  and  have  been  incorporated  in  NASTRAN  from  the  very  beginning  of 
the  code. 

The  basic  displacement  assumption  for  TRIARG  is: 


u  =  ql  +  q2r  +  q3z 


w  =  q4  +  q5r  +  qgz 


(6  terms) 


*»» 


r,u,  Fr 


But  the  generalized  harmonic  displacements  for  TRAPRG  are: 


q,  +  q.  r  +  q,  2  +  q.  rz 
In  2n  3n  4n 


w  =  +  q,  r  +  q„  z  +  q„  rz 

n  5n  6n  7n  8n 


(8  terms) 


Stresses  for  TRAPRG  are  "evaluated  at  the  four  nodal  regions  as  well  as  at  a 
fifth  region  which  corresponds  to  a  point  that  is  obtained  by  averaging  the 
coordinates  of  the  four  nodal  points."  (Currently,  many  other  codes  already 
offer  more  accurate  axisymmetr ic  elements  with  mid-side  nodes.) 

Any  point  within  a  solid-of-revolution  element  can  be  located  by  specifying 
a  radial  coordinate  (r) ,  an  axial  coordinate  (z) ,  and  a  circumferential 
position  (4>)  .  In  TRAPRG  and  TRIARG,  the  point  "i"  displacements  U  ^  and  U  ^ 
are  constants  with  respect  to  $.  In  TRAPAX  and  TRIAAX,  the  displacements 
Ur.  ,  uzi»  and  Upjy  as  well  as  the  forces  F^,  F^,  and  at  each  point  i 
are  assumed  to  be  sinusoidal  functions  of  the  circumferential  location  4>. 

Directly  applied  nonaxisymmetr ic  loads,  pressure  loads,  gravity  loads,  and 
thermal  loads  are  applied  to  the  TRAPAX  and  TRIAXX  elements  with  the  same 
method  as  in  the  conical  shell  element  CONEAX. 

Two  shortcomings  are  noted  in  the  documentation  of  solid-of-revolution 
elements  in  NASTRAN.  First,  there  is  no  explanation  or  justification  for  the 
trapezoidal  requirement,  that  is,  two  opposite  sides  of  TRAPRG  must  be  parallel. 
This  strict  trapezoidal  requirement  exists  only  in  COSMIC/NASTRAN ,  and  is  extremely 
restrictive  on  the  analyst's  choice  of  mesh  to  the  point  of  inconvenience.  Second, 
no  mention  is  made  of  the  relationship  or  similarity  of  the  NASTRAN  solid-of- 
revolution  elements  to  those  in  the  SAAS  code.  The  SAAS  II  and  III  codes  (and 
the  asymmetric  ASAAS  code)  also  utilize  the  Wilson  element,  handle  orthotropic 
behavior,  and  have  been  extensively  used  in  the  aerospace  industry  for  ten  years. 
Moreover,  SAAS  offers  the  user  a  choice  of  plane  strain,  plane  stress,  or 
axisymmetr ic  analyses,  a  choice  not  found  in  NASTRAN.  SAAS  has  an  excellent 
reputation  as  a  versatile  and  cheap  code  to  run.  It  seems  to  the  authors  that 
before  and  during  implementation  of  the  NASTRAN  solid-of-revolution  elements, 
there  should  have  been  several  test  cases  made  to  compare  NASTRAN  efficiency  with 
SAAS  III,  and  to  possibly  incorporate  SAAS  numerical  algorithms  and  its  powerful 
Laplacian  mesh  generation  capability.  The  fact  is  that  an  user  who  has  SAAS 
in-house  will  probably  never  use  NASTRAN  for  axisymmetric  analysis  of  solids. 


7-45 


7.6.2.  Advanced  Evaluation 

The  benchmark  problem  chosen  to  evaluate  TRAPRG  efficiency  and  convergence 
character istics  is  a  thick-walled  cylinder  (plane  strain)  under  internal  pressure. 
This  problem  was  selected  because  its  elasticity  solution  is  well-known,  and  it 
has  been  used  by  past  researchers,  including  Clough  and  Rashid  (1965)  and  Wilson 
(1965) .  The  theoretical  equations  for  radial  and  hoop  stresses  and  for  radial 
displacement  are  (Ref.  47) : 
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Since  there  is  no  input  option  to  apply  pressure  on  the  element  face  for 
TRAPRG  in  COSMIC/NASTRAN,  we  had  to  calculate  and  input  equivalent  nodal 
forces.  TRAPRG  was  evaluated  using  three  meshes:  2,  5,  and  10  elements  through 
the  thickness.  Results  are  compared  with  theory  in  Figures  7.6-1  and  7.6-2. 

The  correlation  is  good. 

TRAPAX  was  also  evaluated  using  the  same  three  meshes,  and  the  excellent 
pressurized  thick-walled  cylinder  results  are  plotted  in  Figures  7.6-3  and  7.6-4. 
Unlike  TRAPRG,  TRAPAX  does  have  an  input  option  for  pressure  directly  applied 
to  an  element  face.  We  were  primarily  interested  in  verifying  the  accuracy  and 
correctness  of  coding  of  TRAPAX.  The  consistency  of  results  between  TRAPRG 
and  TRAPAX  verified  this  coding.  Harmonic  surerposition  used  to  handle  non- 
axisymmetr ic  loads  was  not  evaluated  in  this  soudy,  since  this  procedure  is 
fairly  standard  in  all  the  well-known  codes.  Another  reason  TRAPAX  was  not 
evaluated  for  nonaxisymmetric  loadings  was  the  difficulty  of  formulating  a 
pathological  test  case  for  such  loadings. 

7.6.3.  Conclusions 

For  the  "uniform-dilatation"  thick-walled  cylinder  benchmark,  our  evaluation 
results  showed  both  TRAPRG  and  TRAPAX  to  perform  well  even  for  a  crude  2-element 
model.  The  only  drawbacks  are  the  TRAPRG  trapezoidal  modeling  restriction,  lack 
of  a  direct  pressure  load  input  capability  for  TRAPRG,  and  the  absence  of  one 
element  in  NASTRAN  to  handle  plane  strain,  plane  stress,  or  axisymmetric  analyses. 


THEORY 


TRAPAX  VS. 


t 
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7.7.  The  HEXA1 ,  HEXA2,  IHEX1,  and  I HEX 2  Solid  Elements 

COSMIC/NASTRAN  has  four  constant  strain  solid  elements  of  various  shapes 
(TETRA,  WEDGE,  HEXAl ,  HEXA2)  and  three  isoparametric  solid  elements  of  8,  20, 
and  32  nodes  (IHEX1,  IHEX2,  IHEX3) .  We  decided  to  evaluate  only  four  of  these, 
and  compare  their  accuracy  and  performance  together  in  two  test  problems:  a 
slender  cantilevered  beam  under  four  different  loading  conditions;  and  a  pressurized 
thick-walled  cylinder.  In  addition  to  these  two  static  test  problems,  we  also 
studied  the  eigenvalues  of  a  one- inch  cube,  using  the  available  eigenvalue 
extraction  schemes. 

7.7.1.  Theory 

7. 7. 1.1.  HEXAl  and  HEXA2  Constant  Strain  Hexahedral  Elements 

The  HEXAl  and  HEXA2  solid  elements  are  comprised  of  five  subtetrahedra  and 
ten  overlapping  subtetrahedra  respectively,  each  of  which  is  merely  a  three- 
dimensional  equivalent  of  the  constant  strain  triangle.  Each  of  the  eight 
nodes  in  the  hexahedron  thus  consists  of  three  translational  degrees  of  freedom, 
for  a  total  of  24  DOF's.  These  elements  are  subject  to  the  following  restrictions: 
constant  strain  in  each  subtetrahedron;  uniform  and  isotropic  material  properties ; 
uniform  temperature  in  each  tetrahedral  subelement;  translational  degrees  of 
freedom  only  at  each  node;  and  no  capability  for  differential  stiffness,  buckling, 
and  piecewise  linear  analyses.  The  assumed  displacement  field  at  each  node  is 
therefore  linear  in  each  Cartesian  coordinate: 


U  =  ql  +  q2  +  q3y  +  q4Z 


v  =  q5  +  qgx  +  q?y  +  qgZ 
w  =  q9  t  qlQx  +  quy  +  ^z 


2 


This  assumed  displacement  field  produces  a  constant  (uniform)  strain  and 
stress  in  each  element. 

A  hexahedron  can  be  cut  into  five  subtetrahedr a  into  only  two  different 
ways.  HEXAl  represents  a  single  subdivision  into  five  subtetr ahedra ,  while 
HEXA2  uses  the  average  of  the  results  of  the  two  types  of  subdivisions.  HEXA2 
results  in  symmetrical  deformations  when  symmetrical  loads  are  applied  to  a 
symmetrical  hexahedron,  while  HEXAl  does  not  (Ref.  7) .  This  subtle,  but  im¬ 
portant,  distinction  renders  HEXAl  to  be  definitely  inferior  to  HEXA2  and  the 
IHEXi  elements,  especially  in  problems  where  the  out-of-plane  deflections  (or 
mode  shapes)  are  important.  The  spurious  deflections  are  discussed  further  in 
Subsection  7.7. 2.1. 

In  actual  practice,  numerous  researchers  have  all  concluded  that  far  too 
many  constant  strain  solid  elements  are  required  in  a  three-dimensional  analysis 
to  obtain  acceptable  results  (see  for  instance.  References  48  to  50) .  From 
a  practical  standpoint,  these  constant  strain  solid  elements  are  obsolete  by 
current  1980  standards,  should  be  omitted  in  a  modern  general-purpose  finite 
element  code,  and  have  been  rendered  inferior  by  the  isoparametric  solid  elements 
since  the  late  sixties  and  the  early  seventies. 

7. 7. 1.2.  IHEXI  and  IHEX2  Isoparametric  Hexahedral  Elements 

The  IHEXi  isoparametric  hexahedral  solid  elements  are  effective  for 
analyzing  a  three-dimensional  continuum.  The  material  must  have  isotropic 
temperature-dependent  material  properties.  "Isoparametric"  means  the  element 
deformations  are  represented  with  the  identical  interpolating  functions  used 
to  define  the  geometry.  COSMIC/NASTRAN  offers  the  user  a  choice  of  three  iso¬ 
parametric  hexahedron  solid  elements:  8-node  (IHEXI),  20-node  (IHEX2)  and  32- 
node  (IHEX3) .  These  correspond  to  assumed  displacement  variations  of  linear, 
quadratic,  and  cubic  shapes. 


LINEAR  (8  nodes) 


QUADRATIC  OR  PARABOLIC 
(20  nodes) 


CUBIC  (32  n 
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Clough  (Ref.  48)  and  subsequent  researchers  have  evaluated  various 
three-dimensional  solid  elements,  concluding  that  the  isoparametric  elements 
are  superior  to  other  solid  elements.  The  8-node  and  20-node  isoparametric 
solid  elements  now  appear  in  all  general  purpose  finite  element  codes.  Ex¬ 
perience  has  shown  that  the  8-node  linear  element  is  recommended  for  problems 
with  high  shear  deformations  and  stresses,  while  the  20-node  quadratic  element 
appears  to  work  best  in  problems  with  plate  or  beam  bending  type  deformations. 

Of  these  three  elements,  the  20-node  quadratic  isoparametric  solid  element 
apparently  offers  the  best  combination  of  accuracy,  efficiency,  and  cost,  and 
may  well  be  the  most  often  used  solid  element  today. 

The  integrals  in  the  isoparametric  element  stiffness,  mass,  and  load 
matrices  are  evaluated  by  the  use  of  numerical  integration.  Gaussian  quadrature 
is  used  in  NASTRAN,  and  the  user  has  a  choice  of  a  2  x  2  x  2  Gaussian  integration 
scheme  or  3  x  3  x  3  integration  scheme.  For  most  problems,  the  2x2x2  scheme 
seems  to  work  best  for  the  linear  element  (IHEXl) ,  and  a  3  x  3  x  3  scheme  works 
best  for  the  two  higher  order  elements.  In  our  advanced  evaluation,  we 
evaluated  only  HEXAl,  HEXA2,  IHEXl,  and  IHEX2. 

7.7.2.  Advanced  Evaluation 

The  two  test  cases  used  to  evaluate  HEXAl,  HEXA2,  IHEXl,  and  IHEX2  were 
a  slender  cantilevered  beam  (with  four  loading  conditions)  and  an  internally 
pressurized  thick-walled  cylinder.  The  cantilevered  beam  problem  is  commonly 
used  to  test  the  ability  of  a  solid  element  to  represent  bending  behavior. 

The  thick-walled  cylinder  problem  enables  us  to  assess  the  behavior  of  the 
solid  elements  in  an  enforced  radial  dilatation  mode.  A  third  test  problem, 
eigenvalue  extraction  of  a  unit  cube,  allows  us  to  verify  both  NASTRAN  solid 
element  theory  and  eigenvalue  extraction  efficiency. 


V 


. 1.  Cantilevered  Beam 

The  four  loading  cases  for  the  cantilevered  beam  are: 

Midside  nodes  for  IHEX2  only 
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Theoretical  and  finite  element  deflections  and  stresses  are  plotted  in  Figures 
7.7-1  to  7.7-4.  It  is  obvious  from  these  plots  that  of  the  four  solid  elements 
evaluated  in  this  bending  example,  the  only  element  which  performed  reasonably 
well  is  IHEX2.  IHEX2  consistently  produced  good  correlation  with  theoretical 
deflections  and  stresses  in  all  four  load  cases. 

The  next  best  solid  element  is  IHEX1,  but  it  gave  deflections  and  stresses 
only  30  to  50  percent  of  those  predicted  by  theory.  HEXAl  and  HEXA2  both 
correlated  very  poorly  with  theoretical  deflections.  This  lackluster  performance 
of  IHEXl,  HEXAl,  and  HEXA2  in  a  beam  bending  test  problem  can  probably  be  attributed 
to  the  fact  we  used  only  one  layer  of  elements  through  the  thickness.  This 
leads  us  to  the  conclusion  that,  for  these  three  solid  elements,  the 
analyst  should  use  more  than  one  element  through  the  thickness  in  a  beam  bending 
problem  before  he  can  hope  for  reasonably  accurate  answers.  However,  we  note 
that  for  a  similar  problem,  MSC/NASTRAN * s  6-  and  16-node  PENTA  and  8-  and  20- 
node  HEXA  isoparametric  solid  elements  as  well  as  STARDYNE's  8-node  isoparametric 
solid  element  all  correlated  fairly  well  with  theoretical  displacements  even 
for  a  beam  model  with  one  layer  through  the  thickness. 

Our  results  for  HEXAl  in  the  cantilevered  beam  benchmark  also  indicated 
the  existence  of  spurious  modes  in  the  lateral  x-direction  when  the  primary 
deflection  occurs  in  the  z-direction.  These  lateral  sidesway  modes  for  HEXAl 
were  attributed  to  its  composition  of  five  subtetrahedra  and  the  basic  unsymmet- 
rical  nature  of  this  composition.  The  magnitude  of  this  sidesway  deflection 
was  approximately  five  to  ten  percent  of  the  z-deflection. 

The  effects  of  modeling  a  cantilevered  beam  with  more  than  one  layer  in 
depth  and  in  width,  and  slender  versus  short  deep  beams  have  already  been  dis¬ 
cussed  elsewhere  by  several  researchers  (e.g.  References  48,50,54).  These 
references  cite  the  superiority  of  the  20-node  isoparametric  solid  element  in 
slender  beam  bending  problems  over  8-node  isoparametric  solid  elements  and 
constant  strain  as  well  as  higher  order  hexahedra.  However,  for  short  deep 
cantilever  beams  where  shear  deformation  is  important,  the  8-node  isoparametric 
solid  is  better  than  the  20-  node  isoparametric  solid  and  constant  strain  or 
higher-order  hexahedra. 


BENDING  STRESS  a,  (10  psi)  DEFLECTION  (INCHES) 
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lOOO* 


FIGURE  7.7-1.  RESULTS  OF  CANTILEVERED  BEAM  TEST  OF  SOLID  ELEMENTS 


CASE  1. .  END  LOAD 


BENDING  STRESS,  o  ,  (10*  psi)  DEFLECTION  (inches) 
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LENGTH  (inches) 


FIGURE  7.7-3.  RESULTS  OF  CANTILEVERED  BEAM  TEST  FOR  SOLID  ELEMENTS 


CASE  3.  INTERMEDIATE  LOAD 


BENDING  STRESS,  a  ,  (10  psi)  DEFLECTION  (inches) 
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NASTRAN  does  not  offer  the  analyst  an  option  to  improve  solution  accuracy 
by  the  use  of  "incompatible  bending  modes"  (Ref.  55).  These  so-called  "bubble 
modes"  have  often  been  used  effectively  in  bending  problems?  they  improve 
microscopic  equilibrium  within  the  element,  though  violating  interelement 
compatibility.  The  use  of  such  a  device  has  been  used  very  successfully  to 
improve  the  general  performance  of  the  8-node  isoparametric  solid  element 
(for  instance,  the  ANSYS  STIF  45  element). 

For  the  NASTRAN  user,  therefore,  the  only  solid  element  which  should  be 
used  in  bending  problems  is  IHEX2,  the  20-node  isoparametric  solid  element. 

Use  of  a  20-node  element  implies  a  more  complex  and  time-consuming  input, 
leading  to  a  higher  probability  of  input  error.  IHEX3  (the  32-node  isoparametric 
solid  element)  was  not  evaluated  in  this  study  but  the  same  comments  apply. 

7. 7. 2. 2.  Thick-Walled  Cylinder 

This  test  problem  is  identical  to  that  used  for  TRAPRG  (Subsection  7.6). 

The  problem  measures  the  ability  of  the  solid  elements  to  represent  a  dilatational 
type  of  deformation  (in  the  radial  direction  only).  The  displacement  results  for 
the  four  solid  elements  are  shown  in  Figures  7.7-5  to  7.7-8,  and  the  radial  and 
hoop  stress  correlations  with  theory  are  shown  in  Figures  7.7-9  to  7.7-12.  In 
this  problem,  three  mesh  sizes  were  used:  2,  5,  and  10  elements  through  the 
cylinder  thickness.  This  time,  the  radial  displacement  correlations  are  usually 
12  percent  or  better,  and  the  stress  correlations  with  theory  were  all  remarkably 
good.  Even  though  the  thick-walled  cylinder  test  problem  is  not  as  severe  a 
test  on  the  solid  elements  as  a  cantilevered  beam,  the  outstanding  performance 
of  all  four  NASTRAN  solid  elements  is  still  remarkable,  even  for  the  2-element 
model. 
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7. 7. 2. 3.  Eigenvalue  Extraction 


The  NASTRAN  solid  elements  HEXA1,  HEXA2  and  IHEX1  were  subjected  to  three 
different  eigenvalue  extraction  methods  in  a  third  type  of  element  test.  The 
use  of  an  eigenvalue  and  vector  solution  has  been  used  to  study  the  convergence 
characteristics  of  an  element  (Ref.  51),  though,  to  our  knowledge,  this  procedure 
has  previously  been  applied  only  to  two-dimensional  elements.  It  has  been 
shown  that,  provided  that  the  monotonic  convergence  conditions  are  satisfied, 
the  finite  element  analysis  overestimates  the  stiffness  of  the  system,  i.e., 
the  calculated  displacements  are  smaller  (in  some  norm)  than  the  exact  displace¬ 
ments.  Hence  when  studying  the  effectiveness  of  an  element,  we  need  to  investigate 
its  stiffness  characteristics,  and  an  element  that  is  less  stiff  will  be  more 
effective.  One  way  of  investigating  the  element  stiffness  characteristics  is  to 
represent  the  element  stiffness  matrix  on  the  basis  of  its  eigenvector,  namely, 
solving  the  eigenvector  problem 

(  [  K]  -  XI)  <),  =  0 

The  results  of  this  test  are  important  in  several  respects.  First,  the 
eigenvalue  extraction  methods  available  in  NASTRAN  may  be  compared  to  each  other 
and  to  a  known  solution.  The  eigenvalues  are  expected  to  contain  rigid  body 
modes  (  three  translations  and  three  rotations)  and  several  repeated  roots,  thus 

posing  a  good  test  of  the  eigenvalue  extraction  methods  for  both  accuracy  and 
efficiency.  Secondly,  the  problem  was  solved  as  if  it  were  a  modal  analysis. 

Two  different  approaches  were  used  in  this  regard.  The  first  method  was  to 
require  the  program  to  lump  the  mass  according  to  the  particular  element 
algorithm.  The  second  was  to  give  the  element  a  zero  density  and  place  lumped 
masses  at  each  node.  Finally,  the  mode  shapes  produced  from  the  eigenvalue 
extraction  should  be  distinct  and  recognizable,  i.e.,  shear,  flexure, 
dilatational,  etc.  Failure  to  produce  distinct  mode  shapes  indicates  that 
spurious  displacement  modes  exist  in  the  element  formulation. 

The  eigenvalues  obtained  for  these  elements  are  shown  in  Tables  7.7-1  through 
7.7-3.  The  columns  labeled  "unit  cube”  give  the  results  for  one  solid  element, 
one  inch  on  each  side.  (A  density  of  1.0  slugs.  Young's  modulus  of  10  x  106 
psi  and  Poisson's  ratio  of  .3  were  used  for  material  properties).  The  first 
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UNIT 

CUBE 

ESS3SS9I 

UNEQUAL  &  LUMPED  MASS 

EXTRACTION  METHOD 

GIVENS 

INVERSE 

GIVENS 

GIVENS 

Mode  7 

707 

707 

707 

596 

Number  8 

707 

707 

707 

599 

9 

821 

821 

818 

741 

10 

821 

821 

821 

744 

11 

821 

821 

825 

744 

12 

1102 

1102 

1100 

882 

13 

1102 

1102 

1104 

882 

14 

1148 

1148 

1147 

882 

15 

1148 

1148 

1148 

939 

16 

1148 

1148 

1149 

942 

17 

1184 

1184 

1185 

944 

18 

1184 

1184 

1185 

964 

19 

1184 

1184 

1185 

973 

20 

1323 

1323 

1323 

1139 

21 

2043 

2043 

2043 

1606 

22 

2152 

2152 

2151 

1607 

23 

2152 

2152 

2152 

1609 

24 

2152 

2152 

2153 

1762 

TABLE  7.7-1.  EIGENVALUES  FOR  HEXA1  ELEMENT  (HERTZ) 


i 


TABLE  7.7-2.  EIGENVALUES  FOR  HEXA2  ELEMENT  (HERTZ) 
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EXTRACTION  METHOD 

UNIT  CUBE 

UNEQUAL 

CUBE 

UNEQUAL  & 

MASS  OF 

LUMPED 

NODES 

GIVENS 

* 

INVERSE 

* 

GIVENS 

*  GIVENS  * 

Mode  7 

856 

+ 

856 

+ 

849 

356 

T 

Number  8 

856 

+ 

856 

+ 

863 

364 

T 

9 

1255 

+ 

1255 

+ 

1255 

486 

H 

10 

1448 

+ 

1448 

+ 

1448 

489 

H 

11 

1448 

+ 

1448 

+ 

1449 

489 

H 

12 

1529 

+ 

1529 

+ 

1528 

622 

FX1 

13 

1529 

+ 

1529 

+ 

1529 

622 

FX1 

14 

1647 

+ 

1647 

+ 

1644 

622 

FX1 

15 

1647 

+ 

1647 

+ 

1651 

718 

FX3 

16 

1670 

+ 

1670 

+ 

1671 

876 

E 

17 

1873 

+ 

1873 

+ 

1873 

876 

SH 

18 

2163 

+ 

2163 

+ 

2163 

879 

FX2 

19 

2163 

+ 

2163 

+ 

2163 

882 

SH 

20 

2239 

+ 

2239 

+ 

2239 

882 

FX2 

21 

2634 

+ 

2634 

+ 

2632 

885 

FX2 

22 

2649 

+ 

2649 

+ 

2648 

885 

SH 

23 

2649 

+ 

2649 

+ 

2654 

885 

E 

24 

2920 

+ 

2920 

+ 

2921 

1591 

B 

* 

Produced  mode  shape  of  type 

indicated  by  letter 

(see  Figure  7. 

7-13). 

+Mode  shape  nondistinct. 


TABLE  7.7-3.  EIGENVALUES  FOR  IHEX1  ELEMENT  (HERTZ) 
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six  eigenvectors  were  calculated  in  all  cases  to  be  essentially  zero  frequency 
(rigid  body  inodes)  and  therefore  are  not  reported.  As  many  as  eight  repeated 
roots  were  found  in  some  cases.  It  is  interesting  to  note  that  the  Givens 
method  was  by  far  the  most  efficient  method  for  this  problem.  The  typical  cost 
of  one  run  was  approximately  $15  using  the  Givens  method.  The  Inverse  method, 
by  comparison,  generally  cost  approximately  $45  per  run.  The  Determinant  method 
of  solution  was  attempted  for  all  three  NASTRAN  elements.  In  all  cases  the  time 
limit  for  the  run  was  exceeded  (a  cost  of  approximately  $150)  and  no  solution 
was  obtained.  The  results  obtained  by  the  Givens  and  Inverse  method  were  essen¬ 
tially  identical,  therefore,  the  Givens  method  was  used  for  the  majority  of  the 
runs  made. 

The  eigenvectors  produced  by  using  the  exact  cube  dimensions  were  not  well- 
defined,  i.e.,  the  mode  shapes  did  not  correspond  to  intuitively  expected 
forms,  such  as  extension,  shear,  etc.  This  is  probably  because  of  the  non¬ 
uniqueness  of  the  eigenvectors  formed  from  the  unit  cube  geometry.  In  order 
to  accentuate  the  shapes,  an  "unequal  cube"  geometry  was  run.  The  x  dimension 
was  made  0.99  inch,  -y  was  made  1.01  inches  while  z  was  kept  at  exactly  one. 

This  model  produced  eigenvalues  which  are  essentially  equal  to  those  produced  by 
the  unit  cube  but,  in  the  case  of  HEXA2,  produced  well-defined  eigenvectors. 

The  shapes  are  shown  in  Figures  7.7-13(a)  and  7.7-13(b).  The  mode  shapes  can 
be  described  as  follows: 


MODE  I 
(Type  T) 


MODE  II 
(Type  E) 


MODE  III 
(Type  FX1) 


A  twist-type  mode  where  one  face  of  the  cube  is  rotated  with 
respect  to  the  opposite  face.  This  mode  appears  to  be  the 
lowest  for  all  elements.  There  are  two  repeated  eigenvectors 
with  this  mode  shape.  This  mode  is  designated  by  a  "T"  in 
the  Tables. 

An  extension  mode  in  which  the  cube  elongates  in  one  direction 
and  is  reduced  in  the  other  two  directions.  The  volume  is 
preserved.  There  are  two  repeated  eigenvectors  and  two  distinct 
mode  shapes  of  this  type.  This  mode  is  designated  by  an  "E" 
in  the  Tables. 

This  appears  to  be  a  flexure  mode  in  which  two  opposite  sides  of 
the  cube  are  transformed  into  rectangles  each  having  the  longer 
sides  perpendicular  to  the  other.  Three  modes  of  this  type  were 
observed.  This  mode  is  designated  by  a  FX1  in  the  Tables. 


I 
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MODE  IV 
(Type  SH) 

MODE  V 
(Type  H) 


MODE  VI 
(Type  FX2) 

MODE  VII 


This  is  a  shear  mode  in  which  opposite  sides  remain  parallel 
but  a  shearing  motion  is  exhibited  in  each  of  three  planes. 

Three  such  modes  were  found.  The  designation  used  in  the  Tables 
is  SH. 

This  mode  was  designated  an  hourglass  mode  (H  is  the  symbol  used 
in  the  Tables).  The  shape  is  characterized  by  non-zero  deflec¬ 
tions  only  in  one  direction,  i.e.,  ux  motion  only,  etc.  Opposite 
corners  of  the  cube  have  identical  motion.  When  viewed  looking 
at  one  edge,  the  deformed  shape  looks  somewhat  like  an  hourglass. 
Three  such  modes  were  found. 

A  second  flexural  mode  was  identified  in  which  the  cube  dilated 
on  one  face  and  contracted  on  an  opposite  face.  This  was 
designated  as  a  French  Roof  mode  with  the  Table  designation 
FX2.  Three  repeated  eigenvalues  exhibited  this  mode. 

This  mode  is  a  third  flexural  mode  in  which  two  opposite  faces 
each  form  a  diamond  pattern  with  the  long  axis  of  the  diamond 
diagonals  mutually  perpendicular.  Only  one  eigenvalue/eigenvector 
exhibited  this  mode  shape.  It  is  termed  FX3  for  identification 
in  the  Tables. 


MODE  VIII 
(Type  B) 


The  highest  mode  in  all  cases  appears  to  be  a  dilatational  mode. 
This  was  designated  by  "B"  in  the  tables  (for  breathing)  and 
only  one  such  mode  was  found  for  each  element  considered. 


The  last  column  in  the  Tables  gives  the  results  for  each  element  for  the 
case  where  1/8  of  the  element  mass  is  lumped  at  each  node  point.  Several  con¬ 
clusions  can  be  drawn  from  these  runs. 

First,  the  HEXA1  element  did  not  produce  distinct  mode  shapes  for  any  of  the 
conditions  analyzed.  It  should  be  noted  that  this  element  is  composed  of  5 
tetrahedra  and  has  already  been  shown  to  produce  spurious  deflection  shapes 
(Ref.  Section  7. 7. 1.1).  Additionally,  the  masses  are  lumped  at  the  nodes  on 
the  following  basis.  The  mass  calculation  involves  the  calculation  of  the  total 
mass  of  each  tetrahedron  in  the  element  and  assigning  1/4  of  the  total  mass  to 
each  of  the  four  node  points.  This  method  of  distributing  mass  will  not  produce 
eight  equal  masses  at  the  nodes.  This  procedure  results  in  a  dependency  on 
the  node  numbering  scheme,  i.e.,  the  results  will  be  different  for  this  element 
depending  upon  the  order  of  element  node  numbering,  which  is,  of  course,  un¬ 
desirable.  Thus,  the  element  does  not  produce  symmetrically  shaped  eigenvectors 
as  it  should.  Further,  it  is  noted  that  the  eigenvalues  changed  when  the  masses 
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were  lumped  at  the  nodes  by  the  analyst  rather  than  using  the  element  algorithm. 

This  was  not  true  for  HEXA2,  which  uses  the  same  basic  subelement  (tetrahedra) 
and  indicates  that  there  may  be  some  error  in  the  programming  of  the  element. 

Table  7.7-2  shows  the  results  of  these  tests  for  the  HEXA2  element.  It  can  be 
seen  that  the  eigenvalues  are  almost  identical  for  all  cases  run.  This  was 
expected  since  the  masses  are  lumped  equally  at  the  nodes  (because  of  the  10 
overlapping  tetrahedra)  and  the  symmetry  of  the  resulting  displacement  functions 
within  the  element.  It  is  noted  also  that  the  frequencies  are  lower  than  for 
the  HEXA1  element,  which  indicates  that  HEXA2  is  less  stiff  and  thus  should 
produce  better  results. 

Table  7.7-3  shows  the  results  for  IHEX1,  the  eight-node  isoparametric  hexahedron. 
It  was  found  that  recognizable  mode  shapes  were  not  evident  until  the  lumped 
mass  case  was  run.  Even  though  the  documentation  stated  that  a  lumped  mass  for¬ 
mulation  is  used  as  a  default  option  in  the  program,  the  Programmer's  Manual  indi¬ 
cates  that  no  provision  is  made  for  the  computation  of  a  lumped  mass  matrix 
for  the  isoparametric  hexahedron  elements.  Thus,  when  the  eigenvalue  problem 
for  this  element  is  solved  using  the  modal  analysis  rigid  format  and  using  a 
non-zero  material  density,  the  mass  is  calculated  according  to  the  shape  function 
for  the  element.  Therefore  ,  in  order  to  compare  the  element  stiffness  to 
the  HEXAi  elements,  the  lumped  mass  cases  must  be  used.  On  this  basis,  it  is 
seen  that  the  eigenvalues  are  lower  for  IHEX1  than  either  of  the  HEXAi  elements 
indicating  that  the  element  is  less  stiff.  It  is  also  interesting  to  note 
that  the  mode  shapes  were  not  produced  in  the  same  order  as  for  the  HEXA2  elements. 

As  a  further  check  on  this  test,  a  series  of  runs  were  made  using  a  similar 
element  in  the  ANSYS  program,  the  STIF  45  eight-node  isoparametric  hexahedron 
(the  incompatible  bending  modes  which  are  an  option  in  this  element  were  suppressed) . 
The  results  are  shown  in  Table  7.7-4.  As  can  be  seen,  the  eigenvalues  are  very 
close  to  those  found  for  IHEX1  and  the  distribution  of  mode  shapes  was  similar. 

It  is  noted  that  the  method  of  eigenvalue  extraction  used  was  the  Jacobi  method. 


EXTRACTION  METHOD 

UNEQUAL  CUBE 

★ 

UNEQUAL  CUBE  & 

LUMPED  MASS 

* 

JACOBI  1 

JACOBI  1 

Mode  7 

1068 

T 

356 

T 

Number  8 

1094 

T 

365 

T 

9 

1520 

E 

486 

H 

10 

1521 

SH 

488 

H 

11 

1536 

SH 

490 

H 

12 

1537 

SH 

507 

FXl 

13 

1538 

E 

509 

FXl 

14 

1872 

FXl 

512 

FXl 

15 

1873 

FXl 

592 

FX2 

16 

1873 

FXl 

693 

FX2 

17 

2162 

FX3 

698 

FX2 

18 

2536 

H 

721 

FX3 

19 

2536 

H 

878 

E 

20 

2547 

H 

878 

SH 

21 

2642 

FX2 

883 

SH 

22 

2649 

FX2 

887 

SH 

23 

2655 

FX2 

888 

E 

24 

2757 

B 

1589 

B 

* 

Produced  mode  shape  of  type  indicated  by  letter  (see  Figure  7.7-13). 


TABLE  7.7-4.  EIGENVALUES  FOR  ANSYS  STIF  45 
ELEMENT  (HERTZ) . 
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7.7.3.  Conclusions 

Our  evaluation  results  indicate  that  the  overall  quality,  performance, 
and  efficiency  of  the  COSMIC/NASTRAN  solid  elements  are  fair — but  one  must 
be  very  c&  of ul  in  selecting  the  appropriate  element  and  mesh  for  a  particular 
applicatic  . 

The  performance  of  each  solid  element  in  the  evaluation  benchmarks  may  . 
be  summari.  .d  in  the  following  table: 


BENCHM  RK  PROBLEM 

NASTRAN  SOLID  ELEMENT 

HEXA1 

HEXA2 

IHEX1 

IHEX2 

1. 

Cantilevered  Beam 

(5  elenu  its,  1  layer) 

Poor 

Poor 

Poor 

Good 

2. 

Thick-willed  Cylinder 

(2,5,10  elements) 

Good 

Good 

Good 

Good 

3. 

Unit-cube  Eigenvalues 

(1  element,  consistent 

or  lumped  mass) 

Poor 

Fair 

Good 

(Not  eval- 

uated) 

In  the  eigenvalue  extraction  benchmark,  we  also  found  that  the  Givens  or  tridiagonal 
method  was  better  than  the  inverse  power  method  with  shifts,  and  both  were  superior 
to  the  determinant  method  for  this  problem.  (The  FEER  method  was  not  used  in  this 
evaluation,  but  is  expected  to  be  equal  to  or  better  than  Givens.) 

The  fairly  good  performance  of  the  constant-strain  HEXA2  element  in  our  study 
corroborates  with  the  conclusions  reached  in  a  recent  Air  Force  Aero  Propulsion 
Laboratory/NASA  Lewis  joint  study  on  a  three-dimensional  stress  analysis  of  a 
thermally-cycled  double-edge  wedge  geometry  specimen  similar  to  those  used  in 
aircraft  gas  turbine  engines  (Ref.  52).  Drake  et.  al. ,  found  excellent  correlation 
for  both  displacements  and  stresses  between  a  354-element  NASTRAN  HEXA2  model  and 
a  64-element  model  using  a  Garrett/Airesearch-developed  code  called  IS03DQ  which 
used  12-node  isoparametric  solid  elements. 
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7.8.  Minicomputer  Versions  of  NASTRAN 

Minicomputers  offer  many  advantages  to  the  user  and  have  increasingly 
gained  popularity  among  structural  analysis  code  users  and  developers.  Among 
these  advantages  are:  vastly  lower  cost  per  problem  (sometimes  by  an  order  of 
magnitude  or  more) ;  practical  decentralization  and  greater  flexibility;  suit¬ 
ability  for  interactive  graphics;  economical  time-sharing  by  many  users  (up  to 
as  many  as  63  simultaneous  users) ;  shorter  and  predictable  turnaround  time; 
local  control  over  computing  resource  allocation;  and  for  the  first  time,  a 
chance  for  the  software  developer  to  request  operating  system  characteristics 
on  minicomputers  to  be  tailored  to  his  needs.  Typically,  minicomputers  cost 
about  a  quarter  of  a  million  dollars,  contain  virtual  memories,  and  offer  32- 
bit  architecture.  Users  of  minicomputers  usually  rave  about  the  reduced  cost, 
versatility,  turnaround,  and  efficiency. 

Many  U.  S.  minicomputer  manufacturers  compete  for  the  market.  The  leader 
is  Digital  Equipment  Corporation  (DEC) ,  with  approximately  40%  of  the  market 
(Ref.  4) .  Other  well-known  mini  manufacturers  include:  Hewlett-Packard,  Data 
General,  Honeywell  Information  Systems,  Prime,  Modcomp  and  Perkin-Elmer .  The 
1979  revenues  to  U.  S.  minicomputer  makers  were  estimated  to  be  $4.7  billion, 
and  in  1980  135,000  minis  are  expected  to  be  shipped.  Codes  already  on  minis 
include:  SAP 6  and  SAP7 ,  SDRC/SUPERB  and  SUPERTAB,  Westinghouse/WECAN,  ANSYS, 
PDA/PATRAN,  MARC/PM  NASTRAN  and  MENTAT,  and  of  course,  versions  of  NASTRAN. 

To  the  best  of  our  knowledge,  three  minicomputer  versions  of  NASTRAN 
currently  exist: 


NAME 

CODE  DEVELOPER 

MINI 

MINI  MANUFACTURER 

MSC/NASTRAN 

MacNeal-Schwendler  Corp. 

Los  Angeles,  CA 

DEC  VAX- 11/7 80 

Digital  Equipment 
Marlboro,  Mass. 

UA I /NASTRAN 

Universal  Analytics  Inc. 
Playa  del  Rey,  CA 

Perkin-Elmer 

3240 

Perkin-Elmer  Corp. 
Oceanport,  N.  J. 

PM/NAS TRAN 

MARC  Analysis  Research  Corp. 
Palo  Alto,  CA 

Prime  450,550, 
650,750 

Prime  Computer  Inc. 
Wellesley  Hills, 
Mass . 

The  minicomputer  version  of  MSC/NASTRAN  is  most  commonly  available  on  the  DEC 
VAX-11/780,  and  thus  VAX/NAS  TRAN  is  by  far  the  most  popular  minicomputer  version 
of  NASTRAN  around.  DEC  is  also  the  leader  in  implementing  structural  analysis 
and  design  codes  on  minis.  Its  latest  "Engineering  Systems  Software  Referral 
Catalog"  (5th  edition,  1980)  features  no  less  than  forty  structural  analysis  and 
design  packages  available  on  various  DEC  minicomputers.  UAI/NASTRAN  and 
PM/NASTRAN  are  relative  newcomers  to  the  field.  UAI/NASTRAN  is  being  offered  on 
the  Perkin-Elmer  3240,  the  latest  in  the  Series  3200,  which  reputedly  performed 
faster  than  the  DEC  VAX-11/780  in  44  out  of  58  benchmark  tests.  PM/NASTRAN  is 
offered  by  MARC  Analysis  on  Prime  minicomputers,  and  its  Version  80  features 
four  new  elements:  HEXA3  and  HEXA4  (8-node  constant  strain  and  20-node  linear 
strain  isoparametric  solid  elements) ,  and  SHELL1  and  SHELL2  (8-node  linear  strain 
and  12-node  quadratic  strain  doubly-curved  isoparametric  shell  elements) .  To 
our  knowledge,  there  are  no  minicomputer  versions  of  COSMIC/NASTRAN. 
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7.9.  Preprocessors  and  Postprocessors  for  NASTRAN 

One  weakness  of  COSMIC/NASTRAN  is  its  lack  of  a  built-in  interactive 
mesh  generation  capability.  In  the  past  five  years,  many  software  packages 
have  been  developed  commercially  to  facilitate  preprocessing  and  postprocessing 
for  NASTRAN.  These  usually  offer  the  user  a  capability  to  generate  and  modify 
the  model  in  a  time-sharing  mode,  and  an  instantaneous  check  of  model  connectivity 
on  a  cathode-ray-tube  terminal  screen.  In  addition,  several  U.  S.  government 
agencies  and  aerospace  companies  -  such  as  the  Naval  Ship  Research  and  Development 
Center,  NASA  Goddard  Space  Flight  Center,  and  Rockwell  International  Corporation  - 
have  developed  pre-  and  postprocessing  packages  (available  from  COSMIC)  to  aid 
the  NASTRAN  user  in  his  model  generation  and  results  interpretation. 

Here,  only  some  of  the  more  popular  commercial  NASTRAN  pre-  and  post¬ 
processing  packages  will  be  discussed.  These  codes  are  typically  available  in 
minicomputer  versions,  and  most  are  also  available  through  the  leading  service 
bureaus.  Table  7.9-1  shows  some  well-known  commercial  NASTRAN  pre-and  post¬ 
processors,  their  developers,  marketing  partner,  and  availability  on  service 
bureaus.  The  acknowledged  leader  in  the  field  is  SDRC/SUPERTAB ,  the  first  inter¬ 
active  preprocessing  code  to  emerge.  For  finite  element  modeling  and  NASTRAN  pre¬ 
processing,  SUPERTAB  has  the  most  users  and  benefits  from  a  marketing  associative 
effort  with  Tektronix,  the  leader  among  CRT  terminal  manufacturers.  FASTDRAW  is 
offered  by  McAuto  and  supports  finite  element  modeling  for  NASTRAN  and  other 
codes.  GIFTS  was  developed  at  the  University  of  Arizona  under  Professor  H.  Kamel; 
its  popularity  appears  to  have  declined  recently.  A  versatile  preprocessor  which 
is  rapidly  gaining  in  popularity  is  PDA/PATRAN,  a  powerful  second-generation 
finite  element  modeling  package  which  offers  color  interactive  graphics.  MENTAT 
and  GRAFAX  are  less  well-known  packages  which  are  beginning  to  build  a  user 
base.  GENFEM  is  a  popular  European  finite  element  code  which  reportedly  can 
also  generate  NASTRAN  models.  In  the  European  market,  it  competes  with  ASKA,  ASAS, 
SESAM69,  BERSAFE,  FINEL,  SAMCEF ,  TPS10,  and  PAFEC75  (Ref.  21).  Preprocessors  and 
minicomputers  are  here  to  stay.  The  NASTRAN  user  should  take  advantage  of  them 
because  the  future  will  bring  more  powerful  packages,  increased  availability  of 
minis  and  graphics  terminals,  and  reduced  costs  in  model  generation  and  post¬ 
processing  . 
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7.10.  Program  Efficiency  Study 

One  of  the  stated  objectives  of  this  project  was  to  evaluate  the  program 
efficiency.  This  is  a  formidable  task  for  a  general  purpose  program  such  as 
NASTRAN  because  of  the  complex  nature  of  the  program.  For  example,  the  main¬ 
frame  version  that  is  utilized  (i.e.,  CDC,  IBM,  etc.)  will  certainly  affect 
the  results  of  such  a  study.  In  each  version  there  are  numerous  differences, 
such  as  overlay  structure,  use  of  double  precision  arithmetic,  type  of  loader, 
availability  of  central  memory,  etc. ,  that  can  affect  running  time  and 
efficiency.  Even  if  the  study  is  restricted  to  a  particular  machine,  the 
program  efficiency  is  an  elusive  commodity  to  assess.  Not  only  are  there  a 
wide  variety  of  rigid  format  solution  procedures  to  be  evaluated,  but  within 
the  rigid  formats  there  are  often  several  options  that  may  affect  the  efficiency, 
i.e.,  selection  of  eigenvalue  extraction  method.  The  use  of  DMAP  and  ALTER 
features  will  also  have  an  effect  (usually  adverse)  on  the  efficiency  of 
NASTRAN.  However,  there  are  some  general  observations  that  can  be  made  con¬ 
cerning  program  efficiency,  some  of  which  were  learned  during  the  course  of 
this  project  and  others  obtained  in  the  course  of  discussions  with  users. 

It  is  generally  accepted  that  NASTRAN  was  designed  to  solve  large  problems. 
Indeed,  some  of  the  problems  that  have  been  solved  using  NASTRAN  probably 
could  not  have  been  solved  using  any  other  computer  program.  Clearly, 
program  efficiency  is  not  an  influencing  factor  in  such  cases.  The  automated 
substructuring  capability  in  NASTRAN  may  be  particularly  efficient  for 
large  problems  in  which  parts  of  the  geometry  are  repeated.  This  capability, 
which  was  recently  added,  is  probably  the  most  convenient  substructuring 
method  that  is  currently  available,  and  requires  a  minimum  of  file  handling 
by  the  user.  Also,  it  is  easy  to  use  and  can  result  in  large  savings  in 
computer  cost  and  engineering  time. 

It  is  generally  felt  that  NASTRAN  solves  medium  to  large  shell  and  beam 
structures  as  well  or  better  than  most  general  purpose  programs.  This  is 
probably  one  of  the  strongest  virtues  of  the  code.  Also,  it  is  generally 
concluded  that  the  rigid  format  solutions  in  NASTRAN  have  been  well  written 
and  work  efficiently.  However,  use  of  the  ALTER  or  DMAP  features  can  greatly 
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increase  running  times.  As  a  particular  example,  one  user  reported  that 
NASTRAN  was  being  used  for  seismic  analysis  of  piping  systems.  DMAP  was 
used  to  calculate  modal  participation  factors  and  to  perform  the  modal 
superposition.  Running  times  appeared  to  be  almost  an  order  of  magnitude 
greater  than  would  be  expected  using  programs  that  are  normally  used  in 
such  analyses. 

In  order  to  obtain  a  quantitative  measure  of  the  efficiency  of  the 
program  for  one  particular  case,  a  series  of  runs  were  made  in  which  a  static 
analysis  was  performed  of  a  flat  plate  with  a  central  point  load.  (These 
runs  were  actually  an  extension  of  the  convergence  test  described  in  Section 
7.2.2  for  QDPLT) .  The  results  of  these  runs  are  shown  on  Figure  7.10-1. 

The  computer  used  was  a  CDC  Cyber  175.  The  costs  were  obtained  from  the 
day  file  from  each  run  and  were  computed  using  an  internal  algorithm  that 
is  used  by  the  commercial  data  center  (United  Computing  Systems)  where  the 
runs  were  made. 

The  largest  problem  run  was  just  over  2000  degrees  of  freedom,  which 
could  be  classified  as  a  "moderate"  size  problem.  The  results  of  this  study 
indicated  that  both  the  cost  and  CPU  time  increased  approximately  linearly 
with  the  number  of  degrees  of  freedom  over  the  problem  size  range  studied. 

A  linear  increase  in  cost  versus  problem  size  implies  that  the  code  will  be 
relatively  more  efficient  for  larger  problems.  The  time  required  to  perform 
a  matrix  decomposition  is  a  complicated  function  of  machine  timing  constants 
and  matrix  topological  parameters  (see  Ref.  13).  However,  an  approximate 
formula  for  the  time  required  for  matrix  decomposition  is: 

T  =  «  M  N  W  2 
2  r  ms 

where  M  is  an  experimentally  determined  machine  dependent  parameter  for  a 
multiply-add  loop,  N  is  the  order  of  the  matrix  and  is  the  root  mean 

square  wavefront.  The  test  problem  utilized  a  square  array  of  elements  and 
thus  the  rms  wavefront  increased  (approximately)  linearly  with  the  number 
of  degrees  of  freedom  in  the  problem.  Thus  it  would  be  expected  that  run 
times  would  exhibit  a  quadratic  increase  with  the  total  number  of  degrees 
of  freedom.  The  linear  relationship  between  run  time  and  degrees  of  freedom 
implies  that  the  solution  algorithm  is  more  efficient  for  larger  problems 
than  for  smaller  ones.  For  a  more  complete  discussion  on  estimates  of 
run  times  and  costs,  see  Section  14  of  Ref.  11. 


FIGURE  7.10-1.  EFFECT  OF  PROBLEM  SIZE  ON  COMPLETION  TIME  AND  COST 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  evaluation  study  of  COSMIC/NASTRAN  reveal  that  the  code 
can  be  given  an  overall  fair  rating  in  its  element  library,  capabilities,  and 
efficiency.  Documentation  is  rated  good  and  in  our  opinion  NASTRAN  has  the  most 
complete  documentation  of  any  general  purpose  program.  There  are  some  areas  where 
the  documentation  could  be  improved,  but  in  general  it  is  quite  good.  The 
strengths  of  NASTRAN  are:  the  DMAP  capability;  versatility  of  analysis  capabilities 
as  offered  in  the  20  rigid  formats;  eigenvalue  extraction  options;  plotting; 
restarting;  and  automated  multi-stage  substructuring  analysis.  Its  weaknesses 
are:  lack  of  a  higher  order  2-D  solid-of-revolution  element  to  handle  plane  strain, 
plane  stress,  or  axisymmetric  cases;  user  support  and  lack  of  a  "hot-line" 
service;  dif f icult-to-read  documentation;  lack  of  an  adequate  nonlinear  analysis 
capability;  a  program  architecture  which  is  unsuitable  for  interactive  use;  and 
a  time-consuming  and  ineffective  error  correction  procedure.  In  the  following 
subsections,  an  assessment  of  the  program  features  and  conclusions/recommendations 
will  be  covered  separately. 

8.1.  Program  Assessment 
.Element  Library 


Element 

Performance 

Rating 

Comments 

QDPLT 

Fair 

Twist  moment  with  one  edge  clamped  caused  problems;  passed 
convergence  tests;  failed  patch  test.  Inferior  to  many 
other  available  plate  bending  elements. 

TRSHL 

Poor 

Convergence  performance  in  Scordelis-Lo  cylindrical  roof 
benchmark  inferior  to  other  commercial  shell  elements. 

TRIM6 

Good 

Performed  well  in  cantilevered  beam  benchmark  in  comparison 
to  QDMEM  and  QDMEMl  elements. 

TRPLT1 

Good 

Compared  well  in  simply-supported  rectangular  plate  bending 
benchmark  with  QUADl,  Melosh,  and  Clough  elements. 

CONEAX 

Good 

Thermal  loading  capability  is  suspect.  Otherwise,  overall 
performance  and  versatility  are  good. 

TRAPRG 

Fair 

Good  performance  in  thick-walled  cylinder  benchmark,  but 
restrictive  in  use  due  to  trapezoidal  element  requirement. 
Needs  plane  strain  and  plane  stress  capabilities. 

TRAPAX 

Fair 

Good  performance  in  thick-walled  cylinder  benchmark,  trape¬ 
zoidal  requirement;  can  handle  nonaxisymmetr ic  loads. 

8.1.  (cont'd). 


Element 

Performance 

Rating 

Comments 

HEXA1 

Poor 

Poor  performance  in  cantilevered  beam  benchmark  and  modal 
analysis;  existence  of  spurious  lateral  displacement;  ade¬ 
quate  convergence  in  thick-walled  cylinder  benchmark. 

HEXA2 

Pair 

Poor  rating  in  cantilevered  beam  benchmark;  adequate  per¬ 
formance  in  thick-walled  cylinder  and  modal  analysis. 

IHEXl 

Fair 

Surprisingly  poor  performance  in  cantilevered  beam  and 
modal  analysis;  good  in  thick-walled  cylinder.  May  need 
"incompatible  bending  modes"  to  improve  bending  capability. 

IHEX2 

Good 

Overall  good  performance  in  cantilevered  beam  benchmark, 
thick-walled  cylinder,  and  modal  analysis. 

.Program  Documentation 

The  manuals  can  be  made  more  readable  and  user-oriented  with  simple  examples, 
summary  tables,  application  and  usage  hints,  cost  estimates,  etc.  The  User’s 
Guide  (Ref.  11)  is  a  good  start. 

.Error  Correction  Procedure  and  User  Support 

These  areas  were  the  subject  of  many  user  complaints.  However,  the  nature  of 
COSMIC/NASTRAN  development  history,  implementation,  and  maintenance  (Section  2) 
seems  to  make  these  areas  difficult  to  change  and  improve.  Recently  imple¬ 
mented  user  support  efforts  at  COSMIC  appear  to  be  improving  in  this  respect. 

.Program  Architecture 

The  modular  concept  behind  NASTRAN  program  architecture  is  flexible,  but  is  based 
on  computer  core  availability  in  the  late  sixties  and  is  now  difficult  to  adapt 
to  a  time-sharing  system  and  interactive  use. 

.  DMAP 

The  matrix  handling  capability  and  user  control  offered  in  DMAP  are  very  flexible 
and  unique  among  general  purpose  codes. 

.Constitutive  Library 

Adequate,  but  deficient  for  general  orthotropic  or  anisotropic  analyses. 


.Eigenvalue  Extraction 

COSMIC/NASTRAN ' s  four  eigenvalue  extraction  methods  are  the  most  offered  in  any 
code.  All  seem  to  work  well,  except  the  determinant  method  (see  Subsection  7. 7. 2. 3). 

.Nonlinear  Capabilities 

COSMIC/NASTRAN  offers  limited  nonlinear  analysis  capabilities,  which  are  weak  com¬ 
pared  to  many  other  nonlinear  codes  currently  available. 

.Special  Features 

Special  NASTRAN  features  -  such  as  static  analysis  with  inertia  relief,  fully- 
stressed  design  optimization,  cyclic  symmetry,  component  mode  synthesis,  represen¬ 
tation  of  part  of  a  structure  by  its  vibration  modes,  analyses  of  control  systems, 
structure/fluid  interaction,  aeroelasticity,  and  acoustic  cavities  -  are  unique 
among  current  competing  codes. 

.Automated  Multi-stage  Substructuring  Analysis 

AMSS  is  flexible  and  powerful.  This  highly-developed  substructuring  capability 
is  matched  only  by  ANSYS  among  current  structural  analysis  codes. 

.Minicomputer  Versions  of  COSMIC/NASTRAN 

None,  but  three  other  minicomputer  versions  of  NASTRAN  exist:  MSC /NASTRAN , 

UA I /NASTRAN ,  and  PM/NASTRAN. 

.Pre-  and  Post-Processors  for  NASTRAN 

Many  interactive  graphics  pre-  and  post-processors  have  been  developed  in  the 
past  five  years  for  NASTRAN.  Among  these  commercial  packages  are:  SDRC/ SUPERTAB , 
MCAUTO/FASTDRAW ,  PDA/PATRAN,  GIFTS,  MARC/MENTAT,  and  AOS/GRAFAX. 

■Cost  and  CPU's  versus  Mesh  Size 

Both  are  approximately  linear  with  mesh  size  (Fig.  7.10). 

8.2.  Conclusions  and  Recommendations 

•NASTRAN  program  architecture  is  flexible,  but  has  not  kept  pace  with  recent 
developments  in  computer  hardware  and  operating  system  software.  However, 
major  changes  in  the  architecture  would  be  prohibitively  expensive  and  would 
require  a  major  rewriting  of  the  program.  The  gains  would  probably  not  justify 
the  expense. 
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NASTRAN  is  designed  foe  large  problems.  Small  or  moderate  size  problems  are 
expensive  compared  to  other  general  purpose  codes.  The  lack  of  self-contained 
mesh  generation  capability  almost  necessitates  the  use  of  a  pre-processor  to 
develop  a  finite  element  mesh.  This  further  reduces  the  usefulness  of  the  program 
for  small  to  medium  size  problems  because  of  the  additional  step  in  converting 
the  output  of  the  mesh  generation  program  into  NASTRAN  data,  learning  how  to  use 
the  mesh  generator,  and  the  additional  control  cards.  The  novice  and  even  the 
common  user  are  often  frustrated  by  this  additional  effort. 

No  interactive  capability  exists.  Many  of  the  current  general  purpose  codes  can 
be  executed  in  an  interactive  mode.  It  is  expected  that  this  trend  will  become 
more  and  more  prevalent  in  the  near  future.  The  present  program  architecture  of 
NASTRAN  would  be  quite  awkward  to  convert  to  an  interactive  mode. 

The  documentation  is  bulky  and  not  user  oriented.  It  could  and  should  be  improved. 

The  matrix  handling  package  (DMAP)  is  very  powerful.  Better  documentation  and 
examples  of  its  use  would  probably  greatly  benefit  the  novice  and  common  user. 

To  improve  the  efficiency  of  the  finite  element  library  and  to  reduce  confusion 
to  users  (especially  relatively  new  ones) ,  some  of  the  older  and  less  effective 
elements  (such  as  TRSHL ,  HEXA1,  HEXA2)  should  be  removed  from  the  program  unless 
they  have  unique  capabilities  that  cannot  be  performed  by  other  elements.  Improved 
elements  are  needed  in  some  areas  (as  discussed  previously).  The  myriad  of  13 
separate  membrane  and  bending  elements  in  COSMIC/NASTRAN  — TRMEM,  TRIM6,  QDMEM, 
QDMEMl ,  QDMEM2 ,  TRPLT,  TRPLT1 ,  QDPLT,  TRBSC,  TRIA1,  TRIA2 ,  QUAD1,  QUAD 2  --  should 
be  revised,  consolidated  and  improved  to  offer  the  analyst  fewer  but  better  choices. 
A  brief  synopsis  of  each  element,  its  capability  and  limitations  should  be  included 
in  the  Users  Manual  so  that  reference  to  the  Theoretical  Manual  is  unnecessary 
for  most  analyses.  Alternately,  the  Users  Manual  could  describe  the  "recommended" 
elements  and  the  experienced  user  would  have  access  to  the  remaining  elements  if 
desired. 

It  is  our  opinion  that  the  only  hexahedral  elements  offered  should  be  isoparametric 
in  formulation,  with  clearly  documented  hints  on  the  merits  of  reduced  integration 
and  incompatible  bending  modes  in  various  applications.  The  poor  performance  of 
IHEX1  on  the  cantilever  beam  problem  could  be  greatly  improved  if  incompatible 
bending  mode  formulation  was  added  to  this  element.  Otherwise,  the  20-node 
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hexahedron  IHEX2  should  be  used  if  only  one  element  through  the  thickness  is 
used  when  bending  loads  are  present  {for  example,  to  represent  thick  shells). 

The  additional  complexity  in  modeling  using  a  20-node  hexahedron  versus  an 
8-node  element  results  in  increased  engineering  time  and  greater  chance  of  model¬ 
ing  errors. 

.NASTRAN  is  primarily  a  linear  code.  The  nonlinear  features  are  relatively 
few  and  not  extensive. 

.An  organization  which  uses  or  intends  to  use  the  public  domain  versions  of 
NASTRAN  should  be  prepared  to  commit  to  a  significant  investment  in  staff 
for  prompt  support,  consultation  and  system  (job  control)  help.  Our  findings 
indicate  that  without  this  level  of  support,  most  organizations  are  frus¬ 
trated  by  the  lack  of  these  necessary  functions. 

.Error  correction  in  COSMIC/NASTRAN  is  time-consuming  and  ineffective.  The 
time  lag  between  when  errors  are  reported  and  corrected  is  such  that  most 
analysts  cannot  wait  for  the  response  from  the  NASTRAN  subcontractor.  Many 
organizations  which  use  NASTRAN  correct  errors  as  they  are  uncovered  and  may 
or  may  not  report  these  changes  to  COSMIC  or  its  subcontractor.  This  process 
results  in  many  (mostly  minor)  variations  of  NASTRAN  among  the  user  community 
and  uniform  quality  control  is  virtually  impossible.  It  is  incumbent  upon 
the  user  or  his  organization  to  verify  the  accuracy  of  the  program  at  each 
installation. 

.In  our  opinion,  NASTRAN  is  a  program  that  should  be  maintained  for  public  usage 
While  the  program  is  in  some  respects  relatively  outdated,  we  feel  that  the 
program  is  still  extremely  versatile  and  serves  a  very  useful  purpose.  One  of 
the  strongest  arguments  in  favor  of  this  is  the  fact  that  NASTRAN  has  been 
used  by  such  a  large  number  of  analysts  over  the  years.  Such  extensive  usage 
is  one  of  the  best  methods  of  verifying  a  program.  Also,  this  large  community 
of  users  is  familiar  with  the  program  and  its  capabilities.  This  body  of 
expertise  is  a  very  valuable  asset.  The  high  degree  of  modularity  of  the 
program  will  allow  for  the  addition  of  new  elements  and  new  analysis  procedures 
without  adversely  affecting  the  existing  program.  Verification  of  new  features 
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is  relatively  simple  compared  to  some  other  general  purpose  codes  because  of 
this  modularity. 

.  The  selection  of  a  general  purpose  program  for  use  by  a  company  or  organization 
must  include  the  following  factors:  (1)  Cost,  (2)  Capability  (including 
element  library  and  analysis  methods),  and  (3)  User  Support.  The  cost  of  a 
program  includes  not  only  the  initial  cost  of  the  source  code  (or  the  royalty 
fee  for  commercial  programs)  and  computer  time,  but  the  cost  of  maintaining 
the  program  on  the  user's  system.  For  example,  implementing  new  features  that 
are  needed,  error  correction,  and  other  maintenance  functions  can  be  time- 
consuming  and  expensive.  Secondly,  the  capability  of  the  program  (Item  2) 
should  not  be  static,  but  must  change  as  new  developments  in  element  theory, 
analysis  methods,  etc.  are  developed.  Thirdly,  user  support  is  very  impor¬ 
tant.  If  possible,  the  principal  developers  of  the  program  should  be  available 
for  consultation  when  problems  arise.  It  is  extremely  inefficient  for  even 
an  expert  user  to  find  and  correct  programming  errors,  and  this  effort  probably 
is  impossible  for  most  users.  It  is  our  opinion  that  all  of  these  functions 
should  be  provided  by  a  single  organization,  if  possible.  Such  an  organization 
is  also  needed  to  plan  the  future  developments  and  additions  that  will  be 
necessary  if  NASTRAN  is  to  survive  as  a  viable  program. 
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